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ABSTRACT 

We present optical and IR integrated colours and SBF magnitudes, computed from 
stellar population synthesis models that include emission from the dusty envelopes 
surrounding TP-AGB stars undergoing mass-loss. We explore the effects of varying 
the mass-loss rate by one order of magnitude around the fiducial value, modifying 
accordingly both the stellar parameters and the output spectra of the TP-AGB stars 
plus their dusty envelopes. The models are single burst, and range in age from a few 
Myr to 14 Gyr, and in metallicity between Z = 0.0001 and Z = 0.07; they combine new 
calculations for the evolution of stars in the TP-AGB phase, with star plus envelope 
SEDs produced with the radiative transfer code DUSTY. We compare these models 
to optical and near-IR data of single AGB stars and Magellanic star clusters. This 
comparison validates the current understanding of the role of mass-loss in determining 
stellar parameters and spectra in the TP-AGB. However, neither broad-band colours 
nor SBF measurements in the optical or the near-IR can discern global changes in the 
mass-loss rate of a stellar population. We predict that mid-IR SBF measurements can 
pick out such changes, and actually resolve whether a relation between metallicity and 
mass-loss exists. 

Key words: stars: AGB and post-AGB — stars: carbon — stars: mass-loss — 
Magellanic Clouds — infrared: stars — stars: circumstellar matter — stars: evolution 
— galaxies: evolution — galaxies: stellar content 
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1 INTRODUCTION. 

Stellar mass-loss is inseparable from stellar evolution and 
death. It is fundamental to inject enriched material into 
the interstellar medium, and hence a major driver of the 
chemical enrichment and evolution of galaxies. Whereas the 
energetic photons emitted by massive, young stars destroy 
molecular gas, the dense outflows of evolved stars return to 
the ISM both molecular gas and dust grains where more H2 
can then form. Mass-loss is especially important at the tip of 
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the asymptotic giant branch (AGB) for intermediate-mass 
stars (IM0 < MzAMS ^ QMq). Consequently, mass- loss 
plays a central role also in shaping the AGB and planetary 
nebula luminosity functions. It determines the white dwarf 
mass spectrum and cooling times, and the minimum super- 
nova (SN) progenitor mass. Mass-loss thus influences the 
rate of SNe types II, lb, and Ic, and possibly determines 
the mas s of SNe Type la progenitors and imp acts their fre- 
quency l|Bowen fc Willsonlll99ll : Iwillsonlbood ). 

Thermally-pulsing AGB (TP-AGB) stars are the most 
luminous red stars in intermediate-age stellar populations, 
and accordingly affect their integrated properties. The rel- 
evance of AGB stars contrasts with our inadequate un- 
derstanding of this evolutionary phase. The TP-AGB has 
proven especially hard to model, on account of the thermal 
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pulses they suffer; tfie convective dredge-up of processed 
heavy elements to the stellar surface; and the ejection of 
the s tellar envelope that ends the phase (jBruzual fc CharlotI 
There is an extensive bibliography d ealing with the 



diffic u lties in treating the TP- AGB (e.g. , Renzini fc Voli 
198ll : llberi fc Renzinil Il983l : iMarigo et all Il996l : iMaraston 



199a 120051 ). In particular, the relation between the funda- 



mental parameters (luminosity, L; mass, M; and metallicity, 
Z) of the central star and the mass-loss rate, M, is not quite 
well understood and, therefore, controversial. Of course, it 
does not help that the dust in the ejected envelopes pre- 
cludes the direct observation of the stellar photospheres. 
Data with relatively good spatial resolution of the dusty 
cocoons themselves, on the other hand, have been available 
only in the last decade, with the arrival of the Infrared Space 
Observatory (ISO) and, presently, the Spitzer Space Tele- 
scope. 

Customarily, M has been d escrib ed by the empiri- 
cal Reimers' law (jReimersI [l975l . Il977l '). written as M = 
rjLR/M; R{L, M, Z) is the stellar radius, and ?7 is a fit- 
ting parameter. Modifications to this law, motivated by 
the large scatter of observationally determined mass-loss 
around the Reimers' relat ion, have been propose d by, 
among others, iBaud fc Habin e (1983'). lBloecke3 (|l995l ). and 
iGroenewegen fc de Jond 1 19941 ). Those working on the de- 
tailed modeling of mass-loss at the TP- AGB have countered 
that, in actuality, stellar luminosity first increases at con- 
stant mass until the star reaches a "clifF in the log M vs. 
log L plane; after this point, mass-loss depends much more 
steeply on stellar parameters than stated by empirical rela- 
tions, such that the stellar envelope is lost nearly exponen- 
tially in time at roughly constant luminosity (Bowen 1995; 
Willson 2000, her Fig. 7, and references therein). In this 
view, empirical relations are the result of very strong selec- 
tion effects, and only reflect the fact that mass-loss rates are 
measurable for just a fraction of the stars undergoing mass- 
loss at any given time. Thus, in the log M vs. LR/M plane, 
the so-called cliff and Reimers' relation are almost coinci- 
dent (see Fig. 8 of Willson 2000). Before reaching the cliff, 
stars have low, unmeasurable mass-loss rates; conversely, the 
stage beyond the cliff is short-lived, and most likely stars 
will be highly obscured behind a dusty envelope. A sample 
of stars wi th observable mass-loss rates should include, ac- 
cording to I Willson (,2000 ), mainly objects within one dex in 
M of the cliff. 

The purpose of this paper is to produce population syn- 
thesis models that include different mass-loss rates in the 
TP-AGB. To this end, we combine stellar population evolu- 
tionary models with theoretical spectral energy distributions 
(SEDs) that take into account the radiative transfer in the 
dusty circumstellar envelopes. To begin with, the SEDs are 
calculated on the basis of the mass-loss rates included in 
the evolutionary tracks; we then produce SEDs for the same 
stellar types, but with mass-loss rates one order of magni- 
tude above and below the rates in the tracks, in order to 
explore the whole range where mass-loss is observable in the 
optical. Finally, we confront the resulting theoretical broad- 
band colours and fluctuation magnitudes with optical and 
near-IR observations of AGB stars and Magellanic star clus- 
ters. 



2 STELLAR POPULATION SYNTHESIS 
MODELS AND M. 

In this paper we use a preliminary version of the 
Chariot & Bruzual (2009; CB09 henceforth) simple 
stellar population (SSP) evolutionary synthesis mod- 
els to compute isochrones in the age range from 
a few Myr to 14 Gyr, and metal (helium) con- 
tent Z{Y) = 0.0001(0.26), 0.0004(0.26), 0.001(0.26), 
0.002(0.26), 0.004(0.26), 0.008(0.26), 0.017(0.26), 0.04(0.30), 
and 0.07(0.34). The CB09 models are formally identical 
to the Bruzual & Chariot (2003; BC03 hereafter) models, 
but include several important improvements. Firstly, up to 
15M0 and for the metal and He contents indicated above, 
CB09 use the tracks from the models with updated input 
physics bv Bertelli et al.. t20P8) . For stars more massive than 
15M0 , in the range from 20 to 120Mp), CB09 use the so- 
calle d Padova 1994 tracks (lAlongi et al 1 19931: Bressan et al.l 



199^; iFagotto et ai]|l994aH bl: lGirardi et al.lll996D . 

Secondly, in the CB09 models used in this paper, the 
TP-AGB evolution of low- and interme diate-mass stars is 
follow ed according to the prescription of lMarigo fc Girardil 
(|2007l ). This semi-empirical prescription includes several im- 
portant theoretical improvements over previous calculations, 
and it has been calibrated using carbon star luminosity func- 
tions in the Magellanic Clouds (MC) and TP-AGB lifetimes 
(star counts) in MC clusters (we refer to the paper by Marigo 
fc Girardi for deta ils). The reader should be aware that 
iBerteUi et al] ||2008|) use a different set of TP-AGB tracks, 
also based on the lMarigo fc Girardil (|2007l ) prescription, but 
extrapolated to different chemical compositions of the stellar 
envelope. Thes e sets of TP-AGB tra cks are un-calibrated, as 
pointed out bv lBertelli et al.| (|2008l ). since no attempt was 
made to reproduce the available observations. CB09 will dis- 
cuss the differences introduced in the evolutionary models 
by the use of the calibrated or the un-calibrated TP-AGB 
tracks from these authors. 

Th us, the CB09 models discussed here use the tracks 
in the iBertelli et al.l (f2008i 'l atlas up to the end of the 
A GB phase, and e x tend these tracks with the results 
of iMarigo fc Girardil (120071 ) to cover the TP-AGB phase. 
lBruzualTll2007l) has shown tha t models computed follow- 
ing the Marigo fc Girardil ([20o3) prescription have brighter 
K-hand magnitudes and redder near-IR colours than other 
models, e.g., BC03, that use a semi-empirical treatment of 
the TP-AGB evolution based on an older empirical calibra- 
tion of the lifetime of these stars, and an educated guess of 
the mass associated to TP-AGB stars of a given luminosityQ 

The Marigo fc Girardi prescription, as implemented by 
CB09, accounts for 15 evolutionary stages in the TP-AGB 
(six in the O-rich phase, six in the C-rich phase, and three 
in the superwind phase). By contrast, the BC03 models in- 
clude only one evolutionary stage at each of these phases. 
The signature of TP-AGB stars, i.e., the red colour of the 
integrated SSP around 1 Gyr, becomes more relevant in the 
CB09 models than in previous computations by BC03 and 



^ We refer to the p apers by Marigo &: Gira rdi and BC03 for 
details. In particular, iGirardi fc Mari^ (|2007|) de rive c onstraints 
on TP-AGB lifetimes, while iMarieofcGirar 3 ||2007|) compare 
their theoretical initial-final mass relation with empirical data of 
white dwarfs in open clusters and binary systems. 
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other authors (e.g.. iMaraston et all bood ) . For Z = 0.008, 
the TP-AGB stars in the CB09 models contribute close to 
a factor of two more light in the K-hand than in the BC03 
models. At maximum, the TP-AGB contributes close to 70% 
of the A'-light in the CB09 model for Z = 0.008, but only 
40% in the BC03 model. The peak K-hand luminosity in the 
BC03 model occurs at around 1 Gyr, whereas in the CB09 
model it stays high and close to constant from 0.1 to 1 Gyr. 
The evolutionary rate is such that the total number of TP- 
AGB stars present in the CB09 1 Gyr isochrone is 3.4 times 
larger than the number of these stars in the BC03 models. 
The TP-AGB stars represent 0.012% of the total number of 
stars in this population at this age in the CB09 model, but 
only 0.0036% in the BC03 model. A 10*' Mq cluster contains 
181 TP-AGB stars in the CB09 model, but only 53 of these 
stars in the BC03 model. The TP-AGB stars are about 0.3 
mag brighter at K in the BC03 than in the CB09 isochrone. 
The net effect of all these factors is an increase of roughly 
90% in the contribution of TP-AGB stars to the t o tal K 
flux of the Z = 0.008 SSP at this age. See iBruzuafl (|2007l ') 
for more details. 

The CB09 isochrones provide, at any given age, the 
number of stars at each of 310 positions in the (LBOL,Teff) 
plane. To each star along the is ochrone, an SED is assigned 
from the IWestera et al.l (|2002l ) stellar library for all stel- 
lar phases, except for the C-rich and superwind stages of 
the TP-AGB; for these we use several options, including 
the DUSTY models already mentioned above. For each of 
the nine stellar metallicities considered, we have calculated 
CB09 models for four possible choices of the SED assigned 
to TP-AGB stars in the superwind phase: (1) the SEDs used 
in the BC03 and CB09 models (standard models hereafter) 
and (2) model spectra computed with the code DUSTY for: 
(a) dusty envelopes that result from the fiducial M during 
these evolutionary phases, (6) dusty envelopes from AI one 
order of magnitude above fiducial, and (c) dusty envelopes 
from M one order of magnitude below fiducial. The input 
spectra for the DUSTY code are the same spectra used at 
these phases by CB09, albeit modified according to the up- 
dated stellar parameters if th e mass-loss-rate is M x 10 or 
M/10. In all cases we use the lChabrie^ l|2003l ) IMF. 

We then compute the time evolution of several prop- 
erties of single-burst stellar populations, in particular their 
integrated fiuxes and fiuctuation luminosities in the V, B, 
R, I, J, H, and Ks bands, as well as in the IRAC and MIPS 
wavelengths observed by the Spitzer Space Telescope (Sec- 
tion |4|0 

Besides predicting stellar time evolution on the 



The spectra of TP-AGB stars have not been updated in the 
preliminary version of the standard CB09 model used here with 
respect to those used in BC03. For C-type TP-AGB stars and 
stars in the superwind phase at the end of the TP-AGB evolution, 
BC03 and CB09 use pe riod-averaged spectra based on models by 
ISchuItheis et aL and observations of Galactic stars. For 

O-rich TP-AGB s tars and stars at the tip of the RGB, BC03 
and CB09 use the IWestera et al. atlas. The final version 

of GB09 will include new spectra for TP-AGB stars. 
^ Strictly speaking, the integrated fluxes and fluctuation lumi- 
nosities of a stellar population at a given time are not the product 
of the mass-loss rate at that time (or at any time) , but rather the 
effect of the total mass-loss up to that point. 



Hertzsprung-Russell diagram, the new iMarigo fc GirardU 
models self-consistently estimate pulsation modes and pe- 
riods, changes in the chemical composition of the stellar en- 
velopes and, most relevant for the present work, mass-loss 
rates owing to the pulsating, dust driven winds of O- and G- 
richAGB stars. In the case of O-rich stars, [Marigo fc Girardil 
l|2007l ) derive luminosities from M, R, and pulsation period 
(P); next, they calculate stellar mass-loss rates according to 
the stars' evolutionary slopes d (lo g M)/ d (log L), on the 
basis of the lBowen fc WillsonI lll99l|) dynamical atmo s pheres 
for Miras including dust PI The lBowen fc WillsonI (|l99ll ) 
models anticipate that mass-loss rates of O-rich stars dimin- 
ish with decreasing metallicity, due to a less efficient dust 
production and a smaller photospheric radius at a fixed lu- 
minosity. To determine mass-loss rates of C-rich stars, lumi- 
nosities are calculated first from stellar temperature (Tcff), 
A/, and P. At low luminosities, mass-loss is driven mostly by 
stellar pulsation, whereas radiation pressure on dust grains 
plays a secondary role. However, at a critical luminosity 
that depends on stellar mass, dust-driven superw inds take 
place. For the superwind phase. iMarigo fc Girardil calculate 
M jTeS, L, M, P, C /O) b ased on the pulsating wind models 
bv lWachter et all (|2002h . although including in addition an 
explicit dependence on the G/O ratio. The mass-loss at each 
evolutionary stage is then taken into account to consistently 
determine the stellar parameters in subsequent phases. 

Examples of predict ed mass-loss rates are presented by 
IMarigo fc Girardil (|2003) in their Figures 15 and 16. They 
also note that, for all cases and G/O values: (a) M is ul- 
timately controlled by the changes in L and Teff linked to 
He-shell fiashes, (6) most of the stellar mass-loss occurs dur- 
ing the high-luminosity but quiescent stages that preceed 
thermal pulses, and (c) the superwind regime is achieved 
during fundamental mode pulsation. 

To date, there are few cases where evolutionary syn- 
thesis models include the effects of mass-loss i n ste l- 
lar spectra. Among these, iLancon fc Mouhcind l|2002l ). 
iMouhcine fc Lanconj (|2003l ). and iMarastonI (|2005|) use av- 
eraged obs erved SEDs, whereas iBressan et all (|l998h and 
iPiovan etTa l. ( 2003) have pioneered the use of analytical re- 
lations for the mass-loss rates and wind terminal velocities 
that then allow them to model the spectra of dusty circum- 
stellar envelopes. We have proceeded in a manner similar 
to these latter works. For each stellar type in the TP-AGB, 
SEDs of star plus envelope were produced with the radiative 
transfer code DUSTY; as was mentioned before, SEDs were 
also processed for the same stellar types, but with mass-loss 
rates one order of magnitude above and below the rates in 
the tracks, in order to explore the whole range where mass- 
loss is observable in the optical, according to Willsorr (200C|). 
A change in M, in turn, entails variations in L and hence 
stellar lifetime, R, envelope and core masses, Teff, pulsation 
period, dust-to-gas ratio, dust composition, and C/O ratio. 
A major challenge for this work is that we are aiming at 
extrapolating mass-loss rates at a large range of metallic- 
ities, but in fact a good calibration of all the parameters 
involved in individual TP-AGB stars exists only for Galac- 



* d (log M)/ d (log L) = 1, and M = 5.67 xlQ-^M© yr'^ 
along the "cliff' line; this rate marks the onset of the "superwind" 
phase. 
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tic and Magellanic TP-AGB stars. Two possible routes are 
open: we can either limit our work to these metallicities, 
or venture to make predictions for lower and higher metal- 
licities, with the clear caveat that most relations have not 
been tested in these conditions. We choose the second, and 
str ive to keep as close as p ossible to the procedures used 
bv lMarigo fc Girardil l|2007l ) to produce the original stellar 
tracks. 

We present our results in the subsequent sections, while 
the interested reader can follow our calculations in detail in 
the Appendix. CB09 models (i.e., with fiducial mass-loss) 
are reported in tables 1 through 4. Optical and near-IR 
colours as a function of age, for different chemical composi- 
tions, are presented in tables 1 and 2, respectively without 
and with dusty envelopes; fiuctuation amplitudes are listed 
in tables 3 and 4. 



3 BROAD-BAND COLOURS. 
3.1 Individual AGB stars. 

Fig. [T] displays the theoretical colour^ of individual stars 
along the 0.2 Gyr (dotted line), 0.5 Gyr (solid line), and 
9.5 Gyr (dashed line) isochrones for populations with four 
metallicities {Z = 0.004, black; 0.008, cyan; 0.017, red; and 
0.04, magenta), and our four choices of spectra for stars in 
the TP-AGB. 

As a first test of the models, we compare the theoretical 
near-infrared (near-IR) broad-band colours to the observed 
two-colour diagram, [H — K] vs. \J — H], of ind ividual AGB 
stars in the compilation bv lPiovan et al.l l|2003l) : the sample 
is shown in Fig. (2] Theoretical colours of individual stars 
with fiducial (thick lines) and high (thin lines) mass-loss 
rates are also shown in Fig. O superimposed on the observed 
AGB sample. The metallicity and age symbols are the same 
as in Fig. [T] 

An examination of the bottom left panel of Fig.[T](SEDs 
with M/10) illustrates that a population where TP-AGB 
stars are nearly devoid of dusty envelopes cannot explain 
the colour range of observed AGB stars. As for the stan- 
dard CB09 isochrones (top left panel of Fig. not even 
the young populations with massive AGB stars reach beyond 
[J — H] ~ 3 and [H — Ks] ~ 2. By contrast, there are stars 
in the sample with both [J — H] and [H — Ks] > 4. Models 
with fiducial M and dusty envelopes (top right panel) cover 
these values comfortably, but only at young ages and with 
solar or supersolar metallicities. Models with M x 10 (bot- 
tom right panel) and high metallicities fit the reddest stars 
regardless of age, whereas for lower Z also young ages are 
required. 

Figures [1] and [2] suggest that some of the stars in the 
sample with [H — Ks] ^ 2 could be young (~ 0.2 Gyr 
old), comparatively massive (A/ > 4Mq at the beginning 
of the TP-AGB phase), stars with subsolar metallicity, go- 
ing through a superwind phase with a mass-loss rate of ~ a 
few X 10~^ — 10** Mq yr~^, or roughly one order of magnitude 
above the fiducial rate. Stars could continue losing mass at 
these rates for ~ 10^ years. 



^ We use Ks or K, depending on the observations the models 
will be compared with. The difference is very hardly noticeable. 



For a more recent example, the 2MASS [H — Ks] vs. 
[J — Ks] integrated colour s of individual AGB ca ndidates 
in the sample published bv ISrinivasan et al.1 (|2009 l) are dis- 
played in Figure [3] Different colours are used for O-rich, 
C-rich, and "extreme" (based on their 2MASS and IRAC 
colours) AGB objects. The colour range is smaller than that 
in Figure (2] This can be better a ppreciated in Figu re B) now 
in the same scale as Figure (2] the lSrinivasan et al.l sample is 
shown as a cloud of gray points, with our models superim- 
posed. Ages and metallicities of the models are indicated as 
in Figures[l]and[2l The agreement with the dat a is similar to 
t hat a chieved for the older sample compiled bv lPiovan et al.l 
(120031). 



3.2 Star clusters. 

Fig. Elpresents theoretical two-colour diagrams, [H — Ks] vs. 
[J - H], for SSPs with different metallicities {Z = 0.0004, 
blue; 0.008, cyan; 0.017, red) and, again, our 4 choices of 
mass-loss and spectra for stars in the TP-AGB. The model 
ages go from 100 Myr to 14 Gyr. The first thing that stands 
out is the range covered by the colour values. As opposed 
to the colours of individual stars (see Figures [T] [51 and 3)), 
the integrated colours of SSPs are in general confined to the 
very small range 0.3 < [J - H] < 1.0, 0.1 «: [i? - Ks] < 0.6, 
with [J — H] varying slightly more than [H ~ Ks], no matter 
what mass-loss rate is used for the models. 

Still, it should be useful to compare the theoreti- 
cal integrated broad-band colours to those of stellar clus- 
ters. Our first test set c ompri se s the MC clusters mea- 
sured by iGonzalez et al.l (|2004l ) : iGonzalez-Lopezlira et all 
(■20051). These authors assembled 8 artificial "superclusters" , 
by coadding data of 191 star clusters in bins with simi- 
lar ages and metallici ties, according to classes I — VII in 
the lSearle et al.l (Il980l ) SWB categorization scheme, plus an 
ultra-young (pre-SWB class) superclusterlfl The purpose of 
such procedure is to reduce the stochastic uncertainty pro- 
duced by the inadequate sampling, in sparse clusters, of stars 
evolving through short evolutionary phases, of which the 
AGB is a prime example. "Superclusters" , therefore, should 
be more appropriate t est objects than indiv i dual star clus- 
ters (see, for examp le, ISantos fc Fr(^ll997l : lBruzual|[200^ : 
ICervifio et al.|[2002l ). In fact, if the assumption is made that 
the numbers of stars in different evolutionary stages have a 
Poissonian distribution, then the theoretical relative errors 



of integrated colours scale as 



-1/2 



where Mtot is the total 



mass of the stellar population (|Cervino et al . 2002). In what 
follows, we will use coloured regions to represent expected 
±1(7 stochastic errors. The MC "supercluster" ages, that go 
fro m ^ 6 Myr to > 10 Gyr, are not the originally adopted 
bv lGonzalez et al. (200 4^. Instead, we use now the updated 



_ 3t al.1 (fi _ ^_ 

eloped bv lElson fc Falll (|l985llT98a ): this parameter relates 
the ages of LMC clusters to their UBV colours. To each su- 
percluster we assign the age that corresponds to the "cen- 
tral" S'-type of its constituents; the error is set to span the 



calibration bv iGirardi et al.l (1199^ of the S'-parameter de- 



^ Individual clusters of each SWB class were centred, sky sub- 
tracted, multiplicatively scaled to a common photometric zero- 
point and dereddened before coaddition. SMC clusters were mag- 
nified to place them at the distance modulus of the LMC. 
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Table 1. Colours for CB09 models with fiducial mass- loss but without dusty envelopes 





B — V 


V — R 


V — I 


V — J 


V — H 


V — Ks 


V — K 






Z 


= 0.017 


Y = 0.26 








0.005 


-0.142 


-0.027 


-0.043 


-0.115 


-0.031 


-0.010 


-0.017 


0.006 


-0.063 


0.062 


0.144 


0.339 


0.622 


0.699 


0.681 


0.007 


-0.059 


0.070 


0.175 


0.443 


0.793 


0.900 


0.878 


0.008 


-0.075 


0.062 


0.184 


0.524 


0.950 


1.097 


1.071 


0.009 


-0.133 


0.020 


0.161 


0.604 


1.146 


1.354 


1.324 


0.010 


-0.138 


0.029 


0.201 


0.721 


1.311 


1.531 


1.499 


0.020 


-0.071 


0.096 


0.348 


0.967 


1.606 


1.816 


1.781 


0.030 


-0.021 


0.127 


0.369 


0.931 


1.528 


1.709 


1.677 


0.040 


0.002 


0.126 


0.335 


0.815 


1.353 


1.510 


1.480 


0.050 


0.035 


0.138 


0.339 


0.790 


1.296 


1.437 


1.409 


0.060 


0.051 


0.143 


0.340 


0.777 


1.264 


1.399 


1.372 


0.070 


0.087 


0.165 


0.369 


0.818 


1.291 


1.425 


1.401 


0.080 


0.097 


0.172 


0.381 


0.834 


1.304 


1.437 


1.414 


0.090 


0.117 


0.186 


0.405 


0.869 


1.338 


1.468 


1.446 


0.100 


0.130 


0.196 


0.421 


0.896 


1.368 


1.498 


1.475 


0.200 


0.141 


0.203 


0.498 


1.402 


2.078 


2.345 


2.314 


0.300 


0.171 


0.214 


0.508 


1.302 


2.004 


2.301 


2.304 


0.400 


0.224 


0.234 


0.543 


1.306 


1.986 


2.269 


2.273 


0.500 


0.275 


0.255 


0.587 


1.368 


2.017 


2.262 


2.248 


0.600 


0.311 


0.274 


0.621 


1.412 


2.067 


2.315 


2.302 


0.700 


0.387 


0.309 


0.683 


1.477 


2.142 


2.392 


2.385 


0.800 


0.431 


0.331 


0.716 


1.500 


2.166 


2.415 


2.412 


0.900 


0.477 


0.355 


0.755 


1.542 


2.209 


2.454 


2.451 


1.000 


0.521 


0.378 


0.790 


1.587 


2.251 


2.486 


2.479 


1.500 


0.723 


0.543 


1.178 


2.339 


3.150 


3.432 


3.415 


2.000 


0.707 


0.499 


1.010 


2.035 


2.762 


3.006 


2.981 


3.000 


0.781 


0.537 


1.066 


2.100 


2.831 


3.069 


3.045 


4.000 


0.832 


0.561 


1.099 


2.142 


2.876 


3.109 


3.085 


5.000 


0.854 


0.571 


1.113 


2.159 


2.894 


3.125 


3.102 


6.000 


0.874 


0.581 


1.124 


2.171 


2.907 


3.135 


3.112 


7.000 


0.885 


0.586 


1.129 


2.170 


2.903 


3.127 


3.103 


8.000 


0.901 


0.594 


1.142 


2.194 


2.930 


3.153 


3.128 


9.000 


0.907 


0.597 


1.145 


2.195 


2.928 


3.149 


3.123 


10.000 


0.922 


0.604 


1.159 


2.217 


2.953 


3.174 


3.148 


11.000 


0.929 


0.607 


1.163 


2.219 


2.955 


3.175 


3.149 


12.000 


0.938 


0.611 


1.170 


2.227 


2.963 


3.183 


3.156 


13.000 


0.949 


0.617 


1.179 


2.240 


2.977 


3.197 


3.171 


13.500 


0.955 


0.619 


1.183 


2.247 


2.985 


3.204 


3.178 



(This table is available in its entirety in the online journal, and at CDS in machine-readable format. Values for solar metallicity and 
helium content are shown here for guidance regarding its form and content.) 



5-types of all the members, plus and minus t he rms disper- 
sion 5 (log t) = 0.14 found by Gira rdi et all ([l995;) for the 
log t - S relation. 

Figure [6] shows five two-colour diagrams, comparing 
model SSPs with data of the superclusters, respectively 
[V - I] vs. [H - Ks], [V - I] vs. [J - Ks], [V - Ks] vs. 
[H - Ks], [V - Ks] vs. [J - Ks], and [H - Ks] vs. [J - H]. 
Near-IR data hav e been taken from the Two Micron All Sky 
Survey (2MASS. [Skrutskie et al.lll997l ): I -band data were 
retrieved fro m the Deep Near-Inf rared Southern Sky Sur- 
vey (DENIS: lEpchtein et"aLlll997l '): and V data come from 
different sources in the l i teratu re, mostly from the compi- 
lation by Ivan den Berghl (jl98ll . see Gonzalez-Lopezlira et 
al. 2005). Near-IR colours for the superc l usters were de- 
ri ved for the first time by Gonzalez et all (j2004h . [V - /] 
bvlGonzalez-Lope zlira et al l|2005l ). and \V — K.,] specifically 
for the present work. We have rederived [1^ — /] and the near- 
IR colours, however, to make sure that the individual cluster 
centres are right, that both an image and stellar photometry 
are available for all clusters included in each supercluster, 
and that the background subtraction is optimallj Whereas 



the near-IR colours were measured in the supercluster mo- 
saics, using circular apertures with r — 1' , colours involving 
V {[V~I], [V-J], [V-H], and [V-Ks]) were first obtained 
for single clusters, using the (diver se) diaphragms and V 
magnitudes from the compilation by Ivan den Berghl l|l98lh . 
and then averaged to derive colours for the superclusters. 
In these cases, the quoted errors equal the dispersion of the 
individual colours, divided by (TV — 1)^^^, with A'^ the num- 
ber of objects in each supercluster. Measured colours for all 
superclusters are presented in Table 5, together with their 
ages, metallicities, and photometric masses. Except for su- 
pe rcluster SWB L masses are always lower than estimated 
bv lGonzalez et al.1 (|2004l ): for types III-VII, this is due more 
to our younger (on average, half as old) adopted ages than to 
changes between the BC03 (used by those authors) and the 
CB09 models. While the new assumed ages of types pre- 
SWB, I, and II are older (on average, twice as old), type 
I's estimated mass is the same and type II's, about 30% 
lighter. The only noteworthy case is the pre-SWB superclus- 
ter: its new estimated mass is only 10% of the one derived 



' We have discovered, for example, that NGC 1854 and 
NGC 1855, both reportedly type SWB II, are actually the same 
cluster! Their putative centres are listed to be 6"1 apart; their re- 



spective 5-parameters are, according to lElson fc Falll l ll985l) . 24 
and 22, a difference that provides an independent estimate of the 
uncertainty in the assignment of 5". 
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Table 2. Colours for CB09 models with fiducial mass- loss plus dusty envelopes 



Age (Gyr) B -V V - R V - I V - J V - H V - Ks V - K 



Z = 0.017 Y = 0.26 



0.005 


-0.142 


0.006 


-0.063 


0.007 


-0.059 


0.008 


-0.075 


0.009 


-0.133 


0.010 


-0.138 


0.020 


-0.071 


0.030 


-0.021 


0.040 


0.002 


0.050 


0.035 


0.060 


0.051 


0.070 


0.087 


0.080 


0.097 


0.090 


0.117 


0.100 


0.130 


0.200 


0.143 


0.300 


0. 174 


0.400 


0.228 


0.500 


0.279 


0.600 


0.315 


0.700 


0.392 


0.800 


0.435 


0.900 


0.482 


1.000 


0.526 


1.500 


0.746 


2.000 


0.711 


3.000 


0.783 


4.000 


0.834 


5.000 


0.856 


6.000 


0.874 


7.000 


0.886 


8.000 


0.902 


9.000 


0.908 


10.000 


0.922 


11.000 


0.929 


12.000 


0.938 


13.000 


0.950 


13.500 


0.955 



-0.027 


-0.043 


0.062 


0. 144 


0.070 


0.175 


0.062 


0.184 


0.020 


0.161 


0.029 


0.201 


0.096 


0.348 


0. 127 


0.369 


0.126 


0.335 


0.138 


0.339 


0.143 


0.340 


0.165 


0.369 


0. 172 


0.381 


0.186 


0.405 


0.196 


0.421 


0.206 


0.480 


0.215 


0.495 


0.233 


0.527 


0.252 


0.557 


0.269 


0.586 


0.301 


0.638 


0.320 


0.664 


0.342 


0.698 


0.364 


0.731 


0.520 


1.077 


0.504 


0.987 


0.537 


1.052 


0.561 


1.095 


0.571 


1.110 


0.581 


1.120 


0.586 


1.124 


0.594 


1.137 


0.597 


1.140 


0.604 


1.155 


0.607 


1.160 


0.611 


1.167 


0.617 


1.176 


0.619 


1.181 



-0.115 


-0.031 


0.339 


0.622 


0.443 


0.793 


0.524 


0.950 


0.604 


1.146 


0.721 


1.311 


0.967 


1.606 


0.931 


1.528 


0.815 


1.353 


0.790 


1.296 


0.777 


1.264 


0.818 


1.291 


0.834 


1.304 


0.869 


1.338 


0.896 


1.368 


1.494 


2.324 


1.325 


2.106 


1.282 


2.021 


1.290 


1.983 


1.325 


2.021 


1.369 


2.062 


1.378 


2.063 


1.408 


2.088 


1.451 


2.122 


2.151 


2.999 


2.081 


2.872 


2.144 


2.926 


2.209 


2.998 


2.215 


2.997 


2.207 


2.977 


2.193 


2.951 


2.211 


2.967 


2.206 


2.956 


2.230 


2.983 


2.233 


2.985 


2.242 


2.994 


2.255 


3.008 


2.261 


3.014 



■0.010 


-0.017 


0.699 


0.681 


0.900 


0.878 


1.097 


1.071 


1.354 


1.324 


1.531 


1.499 


1.816 


1.781 


1.709 


1.677 


1.510 


1.480 


1.437 


1.409 


1.399 


1.372 


1.425 


1.401 


1.437 


1.414 


1.468 


1.446 


1.498 


1.475 


2.581 


2.609 


2.395 


2.425 


2.276 


2.297 


2.172 


2.179 


2.213 


2.220 


2.261 


2.270 


2.265 


2.277 


2.287 


2.298 


2.305 


2.310 


3.192 


3.188 


3.054 


3.045 


3.111 


3.094 


3.181 


3.159 


3.186 


3.163 


3.172 


3.148 


3.147 


3.123 


3.164 


3.140 


3.154 


3.129 


3.182 


3.157 


3.186 


3.160 


3.195 


3.169 


3.210 


3.184 


3.217 


3.191 



(This table is available in its entirety in the online journal, and at CDS in machine-readable format. Values for solar metallicity and 
helium content are shown here for guidance regarding its form and content.) 



bv lGonzalez et all (|2004 ). and hence its SBF amplitudes are 
shown now with significantly larger uncertainties. 

In all the panels, the data (solid black circles with er- 
ror bars) are compared to models with a fiducial mass-loss 
rate and 3 different metallicities {Z = 0.0004, blue; 0.008, 
cyan; 0.017, red), that bracket those of the superclusters 
(0.0007 ^ Z ^ = 0.01; Frogel et al. 1990, assuming that 
Zq = 0.017; Cohen 1982). The expected ±la error bars for 
the models, shown as coloured bands, have been calculated 
as in Gonzalez et al. (2004, Appendix), assuming a stellar 
population of 5 x 10"^ Mq. This value is conservative, and 
representative of the MC superclusters; they have between 
~ 10^ and ~ 3 X 10® Mq. The two superclusters where 
the effects of AGB stars should be more important (aged 
160 and 450 Myr, respectively) have both ~ 5 x lO'^ M© 
in stars. Roughly, the models have no problem explaining 
the colours of the superclusters, although we notice that the 
pre-SWB supercluster lies quite far away from the models 
in the bottom left panel of the figure {[H — Ks] vs. [J — H]). 
We hypothesise that th is su percluster might suffer from ad- 
ditional extinction (see iGonz alcz-Lopczlira ct al. 20051) and 
systematic effects. Firstly, in the optical and even in the 
near-IR, there might be a selection against very young clus- 
ters that are the product of the most energetic star forma- 
tion. In fact, ahhough 2MASS data for NGC 2070 (30 Dor) 



were available, they were not used to build the superclus- 
ter type pre-SWB because they showed abundant nebular 
emission and dust extinction. Also, regarding very young 
clusters, the assumption that the addition of many small 
objects is statistically equal to a large one will fail, if none 
of the small clusters is massive enough to produce the most 
massive stars. This problem is not relevant after a few 10^ 
yr, when these stars die and cease to contribute to the clus- 
ter's light. 

Perhaps more useful to assess the models is the compar- 
ison with the data in the age-colour planes. Figure [7] shows 
the data compared to colour versus age for models with a 
fiducial mass-loss rate. The trends of colour with age shown 
by the data are very closely followed by the models. 

It is also instructive to compare our mode ls to the clus- 
ter sarnple already compiled and presented bv lPiovan et al.l 
12003) Fl Once again, two-colour diagrams are shown in Fig. 
[H [V - K] vs. [H - K]; [V - K] vs. [J - K]; and [H - K] 
vs. [J — H]. The data are clusters mostly younger than 1.5 
Gyr, and are displayed as filled circles, with average error 
bars shown for each panel. The expected ilcr error bars for 
the models are depicted once more as coloured regions, ex- 

* The data were extracted by us from Piovan et al.'s paper with 
Dexter dOemleitner et ahllioOll ). 
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Table 3. SBF amplitudes for CB09 models with fiducial mass-loss but without dusty envelopes 



Age (Gyr) 


Mb 




Mr 


Mj 


Mj 














Z ^ 


0.017 Y ^ 


0.26 








0.005 


-7.637 


-7.809 


-8.061 


-8.340 


-8.999 


-9.510 


-9.644 


-9.617 


0.006 


-7.446 


-7.697 


-8. 168 


-8.678 


-9.705 


-10.389 


-10.559 


-10.521 


0.007 


-7.171 


-7.415 


-7.887 


-8.432 


-9.546 


-10.286 


-10.487 


-10.446 


0.008 


-6.885 


-7.052 


-7.457 


-8.033 


-9.295 


-10. 137 


-10.390 


-10.347 


0.009 


-6.427 


-6.470 


-6. 701 


-7.330 


-8.979 


-9.998 


-10.330 


-10.286 


0.010 


-6.096 


-6.092 


-6.352 


-7. 155 


-8.963 


-9.984 


-10.313 


-10.270 


0.020 


-4.936 


-4.972 


-5.263 


-6.051 


-7.692 


-8.698 


-8.999 


-8.954 


0.030 


-4.390 


-4.500 


-4.818 


-5.491 


-6.951 


-7.945 


-8.232 


-8. 186 


0.040 


-3.884 


-4.049 


-4.370 


-4.973 


-6.338 


-7.338 


-7.621 


-7.573 


0.050 


-3.382 


-3.658 


-4.031 


-4.624 


-5.937 


-6.929 


-7.205 


-7. 156 


0.060 


-2.990 


-3.367 


-3. 795 


-4.406 


-5.737 


-6.742 


-7.027 


-6.977 


0.070 


-2.251 


-2.832 


-3.385 


-4.052 


-5.581 


-6.665 


-7.028 


-6.975 


0.080 


-1.895 


-2.547 


-3. 172 


-3.883 


-5.510 


-6.626 


-7.017 


-6.964 


0.090 


-1.629 


-2.341 


-3.019 


-3. 754 


-5.397 


-6.517 


-6.914 


-6.860 


0. 100 


-1.326 


-2.066 


-2. 799 


-3.577 


-5.291 


-6.428 


-6.837 


-6. 782 


0.200 


-0.291 


-0.671 


-1.485 


-3.217 


-6.837 


-8.026 


-8.527 


-8.478 


0.300 


0.083 


-0.254 


-1. 139 


-2.915 


-6. 166 


-7.531 


-8.121 


-8.157 


0.400 


0.508 


0.117 


-0.786 


-2.666 


-5.760 


-7. 132 


-7.724 


-7.766 


0.500 


0.812 


0.339 


-0.595 


-2.622 


-5.698 


-6.981 


-7.496 


-7.492 


0.600 


0.992 


0.435 


-0.555 


-2.597 


-5.600 


-6.887 


-7.406 


-7.410 


0.700 


1.364 


0.619 


-0.483 




-5.359 


-6.664 


-7. 198 


-7.221 


0.800 


1.552 


0.694 


-0.461 


'247^ 


-5. 199 


-6.523 


-7.072 


-7. 109 


0.900 


1.721 


0.744 


-0.447 


-2.419 


-5.072 


-6.405 


-6.967 


-7.012 


1.000 


1.890 


0.805 


-0.405 


-2.337 


-4.968 


-6.300 


-6.869 


-6.912 


1.500 


1.985 


0.421 


-0.890 


-2.793 


-5.305 


-6.491 


-7.032 


-7.072 


2.000 


2.420 


0.991 


-0.160 


-1.837 


-4.897 


-6.073 


-6.568 


-6.535 


3.000 


2.632 


1.171 


0.045 


-1.563 


-4.603 


-5.795 


-6.317 


-6.305 


4.000 


2.893 


1.321 


0.201 


-1.360 


-4.433 


-5.640 


-6.174 


-6.169 


5.000 


2.948 


1.411 


0.307 


-1.218 


-4.285 


-5.499 


-6.034 


-6.032 


6.000 


2.513 


1.413 


0.384 


-1.076 


-4.138 


-5.358 


-5.885 


-5.880 


7.000 


2.745 


1.540 


0.490 


-0.958 


-4.036 


-5.262 


-5.785 


-5.773 


8.000 


2.694 


1.554 


0.519 


-0.912 


-4.025 


-5.248 


-5.765 


-5.744 


9.000 


2.844 


1.635 


0.593 


-0.838 


-3.995 


-5.220 


-5.734 


-5.705 


10.000 


2.786 


1.640 


0.611 


-0.815 


-3.960 


-5.181 


-5.692 


-5.663 


11.000 


2.783 


1.685 


0.664 


-0.767 


-3.917 


-5.139 


-5.650 


-5.620 


12.000 


2.716 


1.703 


0.696 


-0.732 


-3.880 


-5.102 


-5.613 


-5.582 


13.000 


2.818 


1.744 


0.729 


-0.704 


-3.856 


-5.077 


-5.587 


-5.556 


13.500 


2.788 


1.743 


0.734 


-0.695 


-3.846 


-5.066 


-5.576 


-5.545 



(This table is available in its entirety in the online journal, and at CDS in machine-readable format. Values for solar metallicity and 
helium content are shown here for guidance regarding its form and content.) 



cept that now we are assuming a population of 10^ AIq . The 
models have a fiducial mass-loss rate; and metallicities Z = 
0.008 {cyan) and Z = 0.017 (red), that encompass those of 
the clusters. Model ages go from 100 Myr to 1.5 Gyr. 

The match is comparable to that achieved by Piovan et 
al. (2003, their Figure 16) for most colours, and considerably 
better for [V — K] vs. [H — K], even though our model ages 
in the figure stop at 1.5 Gyr (theirs go up to 15 Gyr). 

We have seen so far that our models are able to fit both 
near-IR colours of most single AGB stars, and integrated op- 
tical and near-IR colours of star clusters with different ages 
and metallicities. Integrated colours of star clusters, how- 
ever, do not seem to be able to discriminate between differ- 
ent choices of global mass- loss rates. Next, we will investi- 
gate whether surface brightness fluctuations are potentially 
sensitive to different mass-loss rates in stellar populations. 



4 SURFACE BRIGHTNESS FLUCTUATIONS, 
METALLICITY, AND MASS-LOSS. 

The technique of surface brig htness fluctuation (S BF) mea- 
surements was introduced by iTonrv fc Schneider! ([l988) as 
a way to derive distances to early-type galaxies. The fluc- 
tuation flux (denoted /) is the ratio between the vari- 
ance and the mean of the stellar luminosity function 



(|Tonrv fc Schneidej 1 19881 : iTonrv et al.l Il990l ). scaled by 
(47rci^)~^, where d is the distance. This is expressed as fol- 
lows: 



1 En J? 
47rd2 Eriik 



(1) 



rii and h are, respectively, the number of stars of type i, and 
their luminosity. 

In the case of galaxies, the fluctuation magnitud^l 
m is measured through the spatial fluctuations in their 
surface brightness, and the distance is found by com- 
paring m with empirically calibrated relations that give 
the absolute fluctuation magnitude, M, in a photomet- 
ric band as a functio n of a certain broadband colour, 
in a given range (e.g., IWorth cv 1993a b; P ahrc fc Mould 
!: iTonrv et al]ll997l: lAihar et al.j[l997l: iLiu et alj bOOoT 



Mei et al 
20051: 



20qT 



Mei et al. 



I Jensen et al.l 120031: Cantiello et al 



Gonzalez-Lopezlira et al. 
' 20061 : iMei et all |2007| : 



2005a||b|; 

'Marm-Franch fc Aparicio _ . 

Cantiello et al. 2007a b; Blakcsle e et all |2009| ). In nearby 
stellar clusters, it is possible to obtain m by performing 
the sums in equation [T] over resolved stars (|Aihar fc Tonrvl 



-2.5 log / -I- zero point. 
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Table 4. SBF amplitudes for CB09 models with fiducial mass-loss plus dusty envelopes 



Age (Gyr) 


Mb 


Mv 


Mr 


M, 


Mj 






Mk 








Z = 


0.017 Y = 


0.26 








0.005 


-7.637 


-7.809 


-8.061 


-8.340 


-8.999 


-9.510 


-9.644 


-9.617 


0.006 


-7.446 


-7.697 


-8. 168 


-8.678 


-9.705 


-10.389 


-10.559 


-10.521 


0.007 


-7.171 


-7.415 


-7.887 


-8.432 


-9.546 


-10.286 


-10.487 


-10.446 


0.008 


-6.885 


-7.052 


-7.457 


-8.033 


-9.295 


-10. 137 


-10.390 


-10.347 


0.009 


-6.427 


-6.470 


-6. 701 


-7.330 


-8.979 


-9.998 


-10.330 


-10.286 


0.010 


-6.096 


-6.092 


-6.352 


-7. 155 


-8.963 


-9.984 


-10.313 


-10.270 


0.020 


-4.936 


-4.972 


-5.263 


-6.051 


-7.692 


-8.698 


-8.999 


-8.954 


0.030 


-4.390 


-4.500 


-4.818 


-5.491 


-6.951 


-7.945 


-8.232 


-8. 186 


0.040 


-3.884 


-4.049 


-4.370 


-4.973 


-6.338 


-7.338 


-7.621 


-7.573 


0.050 


-3.382 


-3.658 


-4.031 


-4.624 


-5.937 


-6.929 


-7.205 


-7. 156 


0.060 


-2.990 


-3.367 


-3. 795 


-4.406 


-5.737 


-6.742 


-7.027 


-6.977 


0.070 


-2.251 


-2.832 


-3.385 


-4.052 


-5.581 


-6.665 


-7.028 


-6.975 


0.080 


-1.895 


-2.547 


-3. 172 


-3.883 


-5.510 


-6.626 


-7.017 


-6.964 


0.090 


-1.629 


-2.341 


-3.019 


-3. 754 


-5.397 


-6.517 


-6.914 


-6.860 


0. 100 


-1.326 


-2.066 


-2. 799 


-3.577 


-5.291 


-6.428 


-6.837 


-6. 782 


0.200 


-0.291 


-0.679 


-1.538 


-2.995 


-7.066 


-8.424 


-8.685 


-8. 739 


0.300 


0.084 


-0.265 


-1.170 


-2.763 


-6.233 


-7.724 


-8.144 


-8.223 


0.400 


0.510 


0.093 


-0.805 


-2.477 


-5.646 


-7. 198 


-7.676 


-7.766 


0.500 


0.816 


0.310 


-0.584 


-2.277 


-5.443 


-6.961 


-7.349 


-7.421 


0.600 


0.999 


0.410 


-0.517 


-2.208 


-5.324 


-6.852 


-7.265 


-7.341 


0.700 


1.382 


0.623 


-0.360 


-2.024 


-5.028 


-6.603 


-7.089 


-7.179 


0.800 


1.582 






-1.898 


-4.839 


-6.457 


-6.997 


-7.096 


0.900 


1.765 


0.802 


-0.216 


-1.800 


-4.702 


-6.346 


-6.913 


-7.013 


1.000 


1.949 


0.876 


-0.152 


-1.708 


-4.601 


-6.245 


-6.801 


-6.890 


1.500 


2.091 


0.443 


-0.750 


-2.336 


-4.957 


-6.392 


-6.875 


-6.950 


2.000 


2.426 


0.986 


-0.211 


-1.620 


-5.062 


-6.408 


-6.671 


-6.702 


3.000 


2.635 


1.169 


0.033 


-1.428 


-4.818 


-6.175 


-6.445 


-6.459 


4.000 


2.895 


1.319 


0.197 


-1.326 


-4.841 


-6.219 


-6.460 


-6.454 


5.000 


2.950 


1.410 


0.306 


-1.186 


-4.658 


-6.035 


-6.300 


-6.292 


6.000 


2.513 


1.413 


0.384 


-1.040 


-4.393 


-5.754 


-6.056 


-6.047 


7.000 


2.745 


1.540 


0.491 


-0.920 


-4.204 


-5.548 


-5.882 


-5.870 


8.000 


2.694 


1.554 


0.519 


-0.875 


-4.145 


-5.473 


-5.816 


-5.802 


9.000 


2.844 


1.635 


0.593 


-0.804 


-4.081 


-5.400 


-5.755 


-5.735 


10.000 


2.786 


1.640 


0.611 


-0.784 


-4.057 


-5.369 


-5.730 


-5.708 


11.000 


2.783 


1.685 


0.664 


-0.738 


-4.022 


-5.333 


-5.701 


-5.677 


12.000 


2.716 


1.703 


0.697 


-0.706 


-3.992 


-5.300 


-5.676 


-5.650 


13.000 


2.818 


1.744 


0.729 


-0.681 


-3.970 


-5.274 


-5.656 


-5.630 


13.500 


2.788 


1.743 


0.734 


-0.672 


-3.958 


-5.258 


-5.645 


-5.618 



(This table is available in its entirety in the online journal, and at CDS in machine-readable format. Values for solar metallicity and 
helium content are shown here for guidance regarding its form and content.) 



1994'; 'Gonzalez et al' '2004'; 'Gonzalez-Lop ezlira et al] l2005l : 
Mouhcinc ct al. 2005; Raimondo ct al. 20ol)! 

It is not hard to see that SBFs convey information about 
stellar populations, akin to integrated photometry and spec- 
tra. However, because of their dependence on the square 
of the stellar luminosity, they are especially sensitive to, 
and can provide additional information about the bright- 
est stars at a particular wavelength and at a given evolu- 
tionary phase of a stellar population. Accordingly, it has 
been suggested recently (Cantiello et al. 2003; Raimondo 
et al. 2005, Gonzalez-Lopezhra & Buzzoni 2008) that SBFs 
can be used to study AGB stars in intermediate-age popu- 
lations and, specifically, to investigate their mass-loss rates. 
These works, though, do not explore a possible intrinsic con- 
nection between metallicity and mass-loss, nor consider the 
impact of extinction by dust in the stellar envelope on the 
detectability of mass-losing stars. 

Surprisingly in a way, the relation between metal- 
licity and mass-loss turns out to be controversial, even 
for the dust-driven winds in the TP- AGB. On th e one 
hand, detailed theoretical models (e.g. IWillsonI |200G| ) pre- 
dict that mass-loss shou ld increase with metallicity, and 
iGroenewegen et al.l (|l995h have found, from fits to 8 — 13/im 
spectra, mass-loss rate ratios of 4:3:1 for three O-rich AGB 
stars with similar periods in, respectively, our Galaxy, the 



LMG, and the Small Magellanic Cloud (SMC)0 More re- 
cently, iKalirai et al.l ([2005, 2008) have found evidence of 
a metallicity dependence of the initial-final mass relation- 
ship (between the mass of a white dwarf remnant and 
its main-sequence predecessor) from spectroscopic observa- 
tions of white dwarfs in open clusters. On the other hand, 
iGail fc Scdlmayr (1986) propose that the mass-loss rate is 
proportional to the ratio r/ucxp, where t is the optical depth 
of a dust-driven wind and iicxp is its velocity. Van Loon 
(2000) derives a metallicity-independent mass-loss rate for 
a sample of dust-obscured C and O-rich AGB stars, also in 
the Milky Way, the LMG, and the SMG. Gorrespondingly, 
Ivan LoonI (|2006l ') argues that both t and Wcxp depend on 
the square root of the dust-to-gas ratio, that presum- 
ably is itself linearly proportional to metallicity, such that 
the dependence of M on Z cancels out. Today, however, we 
know that ^ is not constant with metallicity, and that both 
the dust-to-gas ratio and the dust species in the stellar en- 
velopes vary during evolution for a single star and between 



The present-day [Fe/H] ratios fo r the Sun, and B-typ e stars in 
the L MC and the SMC are ~ 3 :2:1 l|Mokiem et al.ll2007^ : accord- 
ing tolL vubimkov et al.l ||2005| ). the Sun and B-type MS stars in 
the solar neighbourhood have the sam e metallicity. The ratios for 
the Sun and F-type stars are ~ 4:2:1 llRussell fc Besselllll989l) . 
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Figure 1. Theoretical two-colour diagrams, \H — Ks\ vs. [ J — of individual stars along the 0.2 Gyr (dotted line), 5.0 Gyr (solid line), 
and 9.5 Gyr [dashed line) isochrones, for populations with different mass-loss rates and metallicities. Top left: standard CB09; bottom 
left: fiducial M/10; top right: fiducial M; bottom right: fiducial AI X 10. Different colours indicate diverse metallicities, i.e., black: Z = 
0.004; cyan: Z = 0.008; red: Z = 0.017; magenta: Z = 0.04. 



different stars with the same initial metalli city (see, e.g., 
iLebzeher eral]|2006l : IPerrarotti fc Gailll2006l ). 

With the aim of addressing the question of the rela- 
tion between Z and M, we compute the time evolution of 
SBF magnitudes of single-burst stellar populations in the 
B, V , R, I, J, H, and Ks bands, in the nine metallicities 
and helium contents mentioned before. The model fluctua- 
tion luminosity L at each wavelength is calculated with the 
following equation (very similar to eq.[T]): 



where the weight Wi is the number of stars of type i per unit 
mass in the population (set as explained in Section rSTjl . and 

h is the lumin osity of stellar t ype 

Following ICervino et all l|2002D . iGonzalez et"al] (|2004D 
demonstrated that the theoretical relative errors of fluctua- 

— 1/2 

tion magnitudes and colours also scale as M^^^ , if a Pois- 
sonian distribution is assumed for the stellar numbers in 
different evolutionary phases. In the rest of this paper, we 
show calculations for stellar populations with 5 x 10^ Mq; 
this number is representative of the MC "superclusters" . 



TiWil^ The corrected stellar weights were, of course, also applied to 

^ ^ TfWili ' ^'^^ calculation of integrated colours used in Section |3] 
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Figure 2. [H — K] vs. [J — H] integrated colours of individual OH/IR and Mira stars from various sources in the literatu re. Blue solid 
circle s: O-rich stars; green open circles: C-stars. The data were originally compiled and properly corrected for extinction bv lPiovan et al.l 
||2003| ). J, H, and K magnitudes brighter than 7 have typical uncertainties of 0.05 mag. Theoretical colours of individual stars are 
superimposed on the sample of observed AGB stars. Thick lines: fiducial mass-loss rate M; thin lines: high mass-loss rate Mx 10. 
Metallicity and age symbols as in Fig. [T] 



5 SBF RESULTS. 
5.1 Model SBFs. 

Figure[9]shows absolute fluctuation magnitudes vs. log (age) 
for standard CB09 models with difTerent metallicities, from 
Z = 0.0004 (~ 1/40 solar) to Z ^ 0.07 4 times solar). 
Coloured regions delimit expected ± 1 cr stochastic errors for 
a stellar population with 5 x 10^ Mq. In the optical bands 
{B to 7), SBFs grow systematically fainter as metallicity in- 
creases. The reason is that the brightest stars (those that 
dominate the SBF signal) will be cooler (redder) at higher 
metallicities. Another consequence of the same phenomenon 
is that errors grow notably larger with metallicity in B and 
V for ages greater than a few Gyr; a higher metallicity trans- 
lates into fewer hot and blue stars, implying larger stochastic 
errors. A very similar trend with metallicity is seen for all 
the mass-loss rates explored in this paper. 



The main difference with metallicity in the near-IR hap- 
pens at very young ages (~ 10^ yr), between the most metal 
poor models {Z = 0.0004), and the rest. Metal poor popula- 
tions will produce SBFs about 2 mag fainter in the near-IR 
at these ages, as a consequence of having fewer red super- 
giants; for the same reason, the stochastic error is largest 
(~ ±1 mag) for these metallicities and ages. Although it 
is strictly true that SBFs increase with metallicity in the 
J H, and Ks bands for populations ~ 100 - 300 Myr old, 
the sensitivity to Z is smaller at 1 — 2/im. This is because, 
once the RGB is born, it is the most important contributor 
to the integrated light of stellar popula tions at these wave- 
lengths (see, e.g.. iGonzalez et al.ll2004l ): the variance (i.e., 
the SBFs) produced by the brightest red stars (both in the 
RGB and AGB) against this already bright background will 
then be less prominent. The onset of the TP- AGB is clearly 
discernible as a peak at about 10* years. 
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Figure 3. 2MASS [H — Ks] vs. [J — H] integrated colours of individual AGB candidates in the sample of lSrinivasan et"al ] ||2009l '). Blue: 
O-rich stars. Red: C-rich stars. Green: "extreme" AGB stars; these are the most luminous AGB stars, losing the most mass. Typical 
photometric errors are 0.05 mag for sources with Kg ^ 11 mag. 



Figure [TD] displays absolute fluctuation magnitudes vs. 
log (age) for models with Z — 0.008 and different mass- 
loss rates. This metallicity {Zq/2.1) is closest to that of 
the youngest half of the MC star clusters whose data are 
presented below. Coloured regions delimit expected ilcr 
stochastic errors for a stel lar population w i th 5 x W^Mq. 
Contrary to the results of iRaimondo et al] (|2005l ). we find 
no strong trend of SBF brightness with mass-loss rate. The 
biggest difference is between the CB09 models without dusty 
envelopes and the models with dusty cocoons, regardless of 
M. This difference is most noticeable in the /-band, where 
models with dust are about 0.5 mag fainter than CB09 stan- 
dard models at all ages after 10* yr. The dusty models ba- 
sically all fall on top of each other. Other metallicities show 
exactly the same behavior, with the exceptions of (1) the 
already pointed out lower brightness and larger dispersion 
of SBFs values for Z = 0.0004 and Z = 0.001 at lO'^ yr, and 



(2) a large scatter for the higher metallicities (Z = 0.017 to 
Z = 0.07) at ages greater than ~ 2 Gyr, in the B, V, and R 
passbands. In the latter case, the dispersion is caused by the 
lower emission from red giants in these bands; it increases 
with metallicity and decreases with wavelength, and goes 
from ~ ±0.5 mag, for Z = 0.017 at V, to ~ ±4 mag for 
Z = 0.07 at B. 

The reason for t his degeneracy is that the selection ef- 
fects highlighted by IWiUsonI (j2000h are intensified by ex- 
tinction, and further exacerbated when the mass-loss rate 
is changed. The probability of detecting the effects of stars 
shedding their envelopes in an exponential fashion, already 
low owing to the short duration of the phase, will decrease 
if the stars are heavily dust-enshrouded. If the mass-loss- 
rate is modified upward, the intrinsic luminosity of stars in 
the superwind stages will increase, but their lifetimes in the 
phase will hence be even shorter and they will be more ob- 
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Figure 4. Theoretical colours of individual stars are superimposed on the sample of observed AGB stars, now shown as gray points. 
Thick lines: fiducial mass- loss rate Af; thin lines: high mass-loss rate Afx 10. Metallicity and age symbols as in Fig. [T] 



scured. The stars losing the most mass sometimes will not 
survive past the first or second superwind stages. Contrari- 
wise, if mass-loss rate is changed downward, stars will last 
longer and be less obscured, but their luminosity will go 
correspondingly down. 



5.2 SBF data. 

SBF magnitu des of MC star c l usters have been determined 
previou sly by iGonzalez et all (|2004 |2005| ) in the near-IR, 
and by iRaimondo et al.l 1 2005h in the optical. As discussed 
earlier, the former built 8 artificial "superclusters" , whereas 
the latter chose to analyse Wide Field Planetary Camera 2 
(WFPC2) Hubble Space Telescope {HST) V and / data of 
a d ozen populous M C globular clusters. 

IConzalez et al.l derived SBFs for the "superclusters" 
within a radius of 1', as prescribed by equation [T] The nu- 
merator was found by summing the square of the fiux of re- 



solved, bright stars, obtained from the 2MASS Point Source 
Catalog; field contamination was minimized by excluding 
from the analysis stars in the range 12.3 < (Ks )o < 14.3 with 
colours (J - Ks)o > 1.2 or (J - Ks)o < 0.4 (|Ferraro et al.1 
1 19951 . Since the sum in the denominator converges slowly, it 
was computed from the total light detected in the images, af- 
ter removal of the emission in an annulus with 2'.0 < r ^ 2.' 5; 
this was assumed to include the contributions from both sky 
and field stars. Absolute fluctuation magnitudes were as- 
signed taking an LM C distance modulus (m — M)o = 18.5 
iFerrarese et al.ll2000l ). 

As it tu r ns ou t, the colour function used by 
IConzalez et al.l (|2004l ) to select cluster stars is adequate 
for older clusters, but is too red for the 2 or 3 youngest 
superclusters. For this reason, we have recalculated near- 
IR fluctuation values for all the superclusters. This time, 
we have statistically remov ed the field populat ion as per 
the procedure described by iMighell et al.l (|l996l ). We com- 
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Figure 5. Theoretical two-colour diagrams, [H — Ks] vs. [J — H], for SSPs with different metallicities and mass-loss rates. Top left: 
standard CB09; bottom left: fiducial M/10; top right: fiducial M; bottom right: fiducial M X 10. Models span an age range between 100 
Myr and 14 Gyr. Blue: Z = 0.0004; cyan: Z = 0.008; red: Z = 0.017. 



pare, for each supercluster, the [J — Ks] versus Ks dia- 
gram of stars within r = 1' of the supercluster centre 
(i.e., the "cluster region"), with that of stars in an annulus 
with 2^0 < r ^ 2'.5 (i.e., the "field"); the former colour- 
magnitude diagram (CMD) presumably includes both clus- 
ter and field stars, whereas we assume that the latter con- 
tains only field stars. For each star in the cluster region 
with mag Ks ± gks and colour [J — Ks] ±o"[j-x3i, we 
count the number of stars in the same CMD with [J — Ks] 
colours within ±MAX(2(T[j_x3] ,0.100) mag and Ks mag 
within ±MAX(2crx3 ,0.200) mag. We caU this number A^sci. 
We also count the number of stars in the field CMD within 
the same AKs by A [J — Ks] bin determined from the clus- 
ter star. We call this number A^'fld. The probability p that 



the star in the cluster region CMD actually belongs to the 
supercluster can be expressed as: 



P 



1 -MIN 



a(iVfld + 1) 
iVsci + 1 



1.0 



(3) 



where a, in this case 0.44, is the ratio of the area of the 
cluster region (tt arcmin ) to the area of the field region 
(2.25 TT arcmin^). Once p is calculated for a given star, it is 
compared to a randomly drawn number ^ p' ^ 1. If p ^ p', 
the star is accepted as a supercluster member; otherwise, it 
is rejected and considered as a field object. 

The numerator of equation [T] was calculated with the 
decontaminated star lists. We have obtained the denomina- 
tor in the same way as before, i.e., from the total hght within 
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Figure 6. Comparison between models and MC "superclusters" ; two-colour diagrams. Top left: \V — I\ vs. [H — Ks\\ top right: \V — I\ 
vs. [J — Ks\; middle left: \V — Ka\ vs. \H — Ka\; middle right: \ V — Ks] vs. [J — Ka\; bott o m lef t: [H — Ks] vs. [J — H]. Filled circles are 
artificial "superclusters" built by iGonzalez et al.l (|2004| . |2005| ) ; iGonzalez-Lopezlira et al] ll2005h . Coloured regions represent SSPs with 
fiducial M and expected itlo" error bars for 5x10^ Mq. Blue: Z = 0.0004; cyan: Z = 0.008; red: Z = 0.017. Supercluster ages go from 6 
Myr to 10 Gyr; model ages span between 3 Myr and 14 Gyr. 



r = l', after subtracting the light in the annulus between 
2'.0 and 2.'5. The new SBF values tend to be fainter than the 
Gonzalez et al. (2004) ones; both sets of values agree within 
1 a, though, with the exception of the SWB I supercluster, 
whose SBF determinations coincide within 3 crri 

An alternative approach for the decontamination, using the 
field CMD as the reference (i.e., "field" stars are removed from 
the cluster regio n, rather than " cluste r" stars kept), has been used 
and described bv lGallart et al.l ||2003| ') . We have tried this method, 
and obtain the same results. We have also tried using the annulus 
with 2.'5 < r sCSfO as the "field". Once again, the results do not 



We have also derived Mj for the eight superclusters 
using DENIS images, and photometry of the point sources 
from the DENIS database^ and the DENIS P oint Source 
Catalogue towards the Magellanic Clouds (jCioni et al.l 
booo). The procedure is the same as for near-IR fluctuation 
magnitudes, except that [/ — J] versus / diagrams have been 

vary significantly, except for the pre-SWB supercluster, although 
even in this case both new derived sets of SBF values lie within 
the (considerably large) errors. 

" Third release of DENIS dat a. The DENIS C onsortium, 2005; 
|http: / / cdsarc.u-strasbg.fr / viz- bin/Cat?B / denis| 
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Figure 7. Comparison between models and MC "superclusters" : colours vs. log (age). Top left: [V — I]; middle left: [V — Kg]; bottom 
left: [J — H]; top right: [H — Ks]; middle right: [J — Kg]. Coloured regions represent SSPs with fiducial M and expected itlcr error bars 
for 5x10^ Mq, coded as in Figure |6] 



employed for the field decontamination. SBF measurements 
for the superclusters are provided in Table 5. 
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Figure 8. Comparison between models and Magellanic star clusters. Top left: [V — vs. [H — K\] top right: \y — K\ vs. \J — K ]\ 
bottom left: [H — K] vs. [J — H]. Filled circles are clusters compiled from the literature and reddening-corrected bv lPiovan et al. 
Coloured regions represent SSPs with fiducial M and expected ±lcr error bars for 10^ Mq. Cyan: Z = 0.008 {cyan); red: Z = 0.017. Most 
of the clusters are younger than 1.5 Gyr; the models range in age between 100 Myr and 1.5 Gyr. Blue triangles and magenta diamonds 
indicate 0.1, 0.3, 0.5, 0.1, 1.1, and 1.5 Gyr, respectively for Z = 0.008 and Z = 0.017; increasing symbol size represents increasing age. 
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Figure 9. Absolute fluctuation magnitudes vs. log (age) for standard CB09 models with different metallicities. Coloured regions delimit 
expected ±la stochastic errors for stellar populations with 5 X 10^ Mq. Black-dotted: Z = 0.0004; green-vertical-hatched: Z = 0.004; 
blue-left-hatched: Z = 0.008; magenta- solid: Z = 0.017; cyan-horizontal-hatched: Z = 0.04; red-right-hatched: Z = 0.07. 
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Figure 10. Absolute fluctuation magnitudes vs. log (age) for Z = 0.008; models with different mass-loss rates. Coloured regions delimit 
expected ilc stochastic errors for stellar populations with 5 X 10^ A^q. Green-vertical-hatched: standard CB09; blue-left-hatched: fiducial 
Af/10; cyan-horizontal-hatched: fiducial M\ red-right-hatched: fiducial MxlO. 



Table 5. Characteristic parameters oiMagellariic superclusters 



Supcrclustcr Liog age (yr) 



?3 
->- 



a b 

Mass (10° Mq) 


[V - I] 


[V - J] 


[V - H] 


[V - Ks] 


[J - Ks] 


[H - Ks] 


[J - Ks] 


Ml 


Mj 






0.07 ± 0.02 


0.34±0.23 


0.89±0.26 


1.32±0.27 


1.64±0.27 


0.47±0.18 


0.34±0.18 


0.81±0.18 


-6.10±a.43 


-5.86±0.64 


-6.76±0.72 


-6.78±0.73 


0.7 ± 0.1 


0.62±0.18 


1.40±0.21 


2.04±0.21 


2.28±0.22 


0.63±0.14 


0.23±0.14 


0.86±0.14 


-5.79±0.26 


-6.77±0.20 


-7.60±0.19 


-7.77±0.19 


0.7 ± 0.1 


Q.54±Q.24 


1.10±0.24 


1.62±0.24 


1.78±0.24 


0.57±0.13 


0.15±0.12 


0.72±0.13 


-4.76±0.35 


-6.59±0.54 


-7.34±0.46 


-7.58±0.41 


0.4 ± 0.1 


0.49±a.l7 


i.oeio.is 


l.,'34±0.18 


1.68±0.18 


0.51±0.11 


0.18±0.11 


0.51±0.11 


-3.a3±a.2a 


-5.97±0.28 


-7.11±0.27 


-7.46±0.29 


0.3 ± 0.0 


0.51±0.35 


1.13±a.33 


1.65±0.33 


1.86±0.32 


0.64±0.16 


0.19±0.15 


O.S3±0.15 


-2.37±a.l9 


-5.43±0.26 


-6.67±0.26 


-7.12±0.25 


0.5 ± 0.1 


1.00±0.23 


1.77±a.26 


2.38±0.26 


2.61±0.26 


0.66±0.16 


0.26±0.16 


0.92±0.16 


-2.30±0.11 


-4.32±0.16 


-5.70±0.19 


-6.46±0.23 


1.8 ± 0.1 


1.02±0.28 


1.78±0.29 


2.39±0.29 


2.58±0.29 


0.61±0.17 


0.17±0.16 


0.78±0.17 


-1.74±0.18 


-3.69±0.20 


-5.04±0.23 


-5.80±0.29 


2.8 ± 0.3 


1.09±0.33 


1.74±0.36 


2.26±0.36 


2.35±0.35 


0.52±0.22 


0.08±0.21 


0.60±0.22 


-2.65±0.35 


-2.71±0.26 


-3.91±0.36 


-4.48±0.51 



pre . 

I . . . 

II . . 



III . 

IV . , 

V . . 

VI . 
VII. 



6.78±0.62 
7.51±0.32 
7.88±0.25 
8.21±0.29 
8.65±0.36 
9.09±0.29 
9.45±0.28 
9.82±0.29 



o.oio±o.q6^B 

O.OlOiO.Og'^ 
O.OlOiO.oK'^ 
O.OlOiO.O^'' 
3c-3±2c-^ 
4c-3±2c-P 
2c-3±lc-04 
7c-4±4c-H 



" Prom the calibration of the 5-parameter bv lGirardi et al 
Masses from near-IR mass-to-light ratios, CB09 models; errors are equal to the dispersion of the results at J, H, and Ks 
CohenI lll982h 

Froeel et alT l ll990l) . assuming Zq = 0.017. 
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Given the si Knificantly greater depth of their HST data, 
iRaimondo et al.l were able to measure both the numerator 
and denominator of equation[T]using resolved stars. They de- 
termined absolute V and / fluctuation magnitudes with the 
same LMC distance modulus (m — M)o = 18.5. The com- 
parison of the described cluster data with the results of the 
models presented in this work is shown in Figures [TTI to [T4l 

Figures [TTI and [rfFI present V, I, J, H, and Ks SBF 
absolute magnitudes vs. log (age) of young and intermediate- 
age MC clusters. Optical {My and A/j) me asurements of in- 
dividu al, rich clusters have been taken from lRaimondo et al.l 
l|2005h . These data are shown as solid circles. For this work, 
we have recalculated near-IR fluctuation values for the su- 
perclusters as described above, and calculated Mi, as well. 
These data are shown as empty triangles. As seen before, 
SBF measurements are not capable of discriminating be- 
tween different mass-loss rates. 

In the near-IR, there is a good overall match between 
models and data. In the optical V and /-bands, there is a 
te ndency for the da t a of s ingle MC clusters first published 
by IRaimondo et al.l (|2005l ) to fall below the model values. 
Among these clusters, the problem is most acute for NGC 
1805 (log t = 7.00 ± 0.20), NGC 1818 (log t = 7 AO ± 0.30), 
and Hodge 14 (log t = 9.30 ± 0.10) In the case of H odge 14, 
severe field contamination forced IRaimondo et al.l to anal- 
yse only a small area (i.e., mass) of the cluster; the resulting 
larger stochastic uncertainty is reflected by the SBF magni- 
tudes error bars. For NGC 1805 and NGC 1818, they men- 
tion a possible 10% incompleteness of the brightest 3 mag 
of cluster stars within 7'.' 2 from the centre, that in principle 
is accounted for by th e shown error ba r s. A l ogical conclu- 
sion -also reached by IRaimondo et al.l (120051 )- is that the 
number of stars in these individual globular clusters is not 
enough to adequately sample the brightest, rarest, TP-AGB 
stars. Additional systematic effects that might impact the 
data of both N GC 1805 and the pre-SWB supercluster of 
IConzalez et al.l have been already discussed in section 15.21 

Results for clusters with Z ~ 0.002 are shown in Fig. 1131 
Figure [T4| compares SBF measurements of old MC clusters, 
with Z ~ 0.0007, to models with Z = 0.001. The agreement 
between models and data is roughly the same if models with 
Z = 0.0004 are used instead. For the one "supercluster" 
shown and contrariwise to what we had seen so far, we find 
that the fit between models and dat a in the optical is better 
than the match in the near-IR. Since IRaimondo et all (|2005l ) 
do not use our same mapping between SWB type and age 
(and thus metallicity) , we show their four oldest clusters 
in both figures [13] and [H Their quoted SWB type is VI, 
like our one supercluster displayed in Fig. [T3] as an empty 
triangle, whereas their listed ages are all over 10 Gyr old, 
i.e., older than our type VII supercluster, and hence we plot 
them too in Fig. (T?) 



6 MID-INFRARED SBF MEASUREMENTS 
AND FUTURE WORK. 

So, can SBF measurements at all, with their sensitivity to 
the brightest stars of a population, provide some insight 
about the mass-loss parameters of unresolved stellar pop- 
ulations? Figure [T^ shows absolute fiuctuation magnitudes 
vs. log (age), again for Z = 0.008, but this time in the mid-IR 
bands observed by the Spitzer Space Telescope. According 
to this figure, in the mid-IR one could begin to distinguish 
intermediate age stellar populations with different mass-loss 
rates, to the point that it might be worthwhile to start ex- 
ploring the effects of other dust mixtures and even of the 
dust chemistry in the stellar envelopes on the integrated 
properties of stellar populations. In a forthcoming paper, 
we will compare models in the mid-IR with Spitzer observa- 
tions of stars in the Large Magellanic Cloud. 



7 SUMMARY AND CONCLUSIONS. 

We have presented optical and IR broad-band colours and 
SBF magnitudes computed from single stellar population 
synthesis models, where the mass-loss rates in the CB09 
evolutionary tracks have been used to produce spectra of 
TP-AGB stars, taking into account the radiative transfer 
in their dusty circumstellar envelopes. Star plus envelope 
SEDs have been processed also for M one order of magnitude 
above and below the fiducial mass-loss rates; for mass-loss 
rates different from the original ones in the tracks, stellar 
configurations and lifetimes have been adjusted. Next, we 
have compared our models to optical and near-IR data of 
single AGB stars and Magellanic star clusters. 

Even though mass-loss regulates lifetimes, luminosities, 
and effective temperatures of stars in the AGB phase, and 
hence their numbers and colours, it turns out that broad- 
band optical and near-IR colours and SBF measurements of 
stellar populations cannot discern global variations in mass- 
loss rate. Worse even than for single stars, the selection ef- 
fects that preclude detection of single stars away from the 
"clifF in a core mass vs. luminosity plot (stars either los- 
ing mass too slowly for the mass-loss to be observed or too 
fast for the star to be detected, due both to lifetime and 
obscuration) make differences in mass-loss rate across whole 
populations completely unnoticeable 

We predict that SBF measurements in the mid-IR could 
begin to pick out intermediate age stellar populations where 
the stars lose mass with rates away from the cliff line, and 
help determine if such deviations actually correlate with 
metallicity. 
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Figure 11. SBFs of Magellanic star clusters I. Observations of young and intermediate-age MC star clusters are compared to models 
with Z = 0.008. ztlcr stochastic uncertainties are shown for s tellar population s with 5 xlO^ Mq, coded as in Figure [TOl Solid circles are 
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Figure 12. SBFs of Magellanic star clusters II. Observations of intermediate-age MC clusters are compared to models with Z = 0.004. 
±1(T stochastic uncertainties are shown for stellar populations with 5 X 10^ coded as in Figure [TOl 



APPENDIX A: MODELLING THE STELLAR 
PLUS ENVELOPE EMISSION. 

Al Mass-loss rate and stellar parameters. 

We treat M as th e independent parameter , and modify L, 
using Figur e 2 in Bow en & Willson (jllU) (reproduced as 
Figure 14 in lMarigo fc G irardi 2007). The former gives lines 
of constant M/M in the log M vs. log L plane, for dusty 
Miras with solar metallicity. 

N ext, w e find the modified stellar radius from the 

llbeni (|1984| ) radius-luminosity-mass relation for evolving 



AGB stars, reproduced as equation 3 in iBowen fc WiUsonl 
|l993)0 and calculate the effective temperature using L — 
B?{T,a/^nQf. 

This is the inve rse of the procedure used by 
iMarigo fc Girardil (|2007l ) to determine M for an oxygen-rich 



16 iJ = 312(L/10'*)0-68{M/1.175)-0-3iS{Z/0.001)0-088('/«p)~''''', 
where L and M are in solar units, (I /Hp) is the ratio of mixing 
length to pressure scale height, and S = for M ^ 1.175 and 
5 = 1 otherwise. 
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Figure 13. SBFs of Magellanic star clusters III. Observations of old MC clusters are compared to models with Z = 0.002. ztlcr stochastic 
uncertainties are shown for stellar populations with 5 X 10^ A^q, coded as in Figure [TOl According to their reported SWB type (VI), 
Raimondo's four oldest MC clusters belong here. 



AGB star (from mass and Z they calculate period, next ra- 
dius, then luminosity, and ultimately mass-loss rate). 

Once the new L, R, and TbS are known, and assuming 
the energy available to a star of type i during an evolution- 
ary phase is constant (i .e., the fuel consumption theorem, 
iRenzini fc BuzzonilllQSd ). the length t2 of the extrapolated 
mass-loss phase can be derived from the relation Liti — 
L2t2, where subscript 1 denotes original parameters and sub- 
script 2 represents the extrapolated ones. If, furthermore, 
we assume that the weight Wi of a stellar type i (i.e., the 
number of stars of initial mass m going through an evolu- 



tionary phase ph, per unit mass of the population) is pro- 
portional to the time spent by such stars in the phase, then 
Wi^2/wi^\ = t2/tl = Ll/I/2. uii,i = 1, by definition, and Wi^2 
is the same for all stars with the same metallicity, once the 
extrapolation factor is fixed, with the following exceptions. 
There are occasions when the product of the mass-loss rate 
and the phase length derived with the above procedure is 
greater than the stellar core mass at the beginning of the 
stage, Mevat. In these cases, Wi^2 = Mevst/M/fa, and the 
weight of the subsequent stage will be zero (i.e., the star will 
not reach the following stage). In other instances, the dust- 
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Figure 14. SBFs of Magellanic star clusters IV. Observations of old MC star clusters are compared to models with Z = 0.001. ztltr 
stochastic uncertainties arc shown for stellar populations with 5 X 10^ A^q, coded as in Figure [TOl According to their reported ages, 
Raimondo's four oldest MC clusters belong in this figure. 



to-mass ratio <& diverges (see ea. lA17l in Section rA2.2l below) . 
consequently the shell opacity also diverges (see ea. lA3l Sec- 
tion |^2]), and hence we assign Wi^2 = to such stars in view 
of their very heavy obscuration. In general, the mass lost in 
any given stage with modified mass-loss rate will be difi'erent 
from the mass lost in the original track, so an adjustment 
to the parameters of the subsequent stage will be needed. 
Logically, we first rectify the stellar mass Mevat and next 
the luminosity L\, moving along a line of constant M /M 
in the aforementioned Figure 2 ofH owen fc WillsonI |l99ll l. 



From then on, we repeat the procedure outlined here: we 
first modify M, and then derive L2, R2, T'e//,2, ^2, Wi^2- 

Fundamental mode (FM) pulsation periods Po also 
change with the stellar configuration . To adjust these, we 
use eq. 12 in lMarigo fc Girardil l|200it ): 



log (Po/days) = - 2.07 + 1.94 log {R/Rq) 
- 0.9 log (M/Mq) (if M < I.5M0) 
= - 2.59 + 2.2 log (R/Rq) 
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Figure 15. Mid-IR absolute fluctuation magnitudes vs. log (age) for Z = 0.008; models with different mass-loss rates. Coloured regions, 
coded as in Figure [TOl delimit expected ztlo" stochastic errors for stellar populations with 5 X 10^ Mq. 



0.83 log (M/Mq) - 0.08 log (Z/IQ-^) 
+ 0.25(F-0.3) (if A/ > 2.5M0). 



(Al) 



Z and Y are, respectively, metallicity and helium content by 
mass; for M between 1.5 and 2.5 Mq, log Pq is interpolated 
linearly, with log M as the independent variable. In the case 
of C-rich stars in the superwind phase, another parameter 
th at changes with M i s the C/0 ratio. For this, we use eq. 23 

■ 'm 



in lMarigo fc Girardil (|2007l )f 



This expression is similar to that proposed bv lWachter et all 



log [A'f/(M0yr~^)] = - 4.529 
6.849 log (reff/2600 K) + 1.527 log (L/W^Lq) 
- 1.997 log (M/Mq) + 0.995 log (Po/650 days) 



||2002|) for winds driven by radiation pressure on dust grains, ex- 
cept that it contains an exphcit d ependence on the C/O ratio. 
Eq. 23 of iMarigo fc Girardil l l2007h has an error: the minus sign 
in front of log Pp sh ould be a plus, as it appears here and in 
IWachter et al.l ||2002|') . 
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+ 0.672 log ■ (A2) 



Eq.miis valid for M ^ 10"^ Mq yr'^ 



A2 Parameters of the dusty envelopes. 

The effect of dust on the stellar SEDs has to be in- 
cluded, both for stars in the original tracks and for those 
whose mass-loss rates and stellar parameters have been 
mo dified. To th i s end , we follow the pr ocedure outlined 
bv lPiovan et al.l hooj) and iMarigo et al.1 (2008) based on 
respec tively, llvezic fc Elitzuij 1 19971 ) and lElitzur fc Ivezid 
l|200ll ). One very important input parameter to produce the 
SED through radiative transfer calculations is the optical 
depth, T, of the envelope. 

For a star of luminosity L and effective temperature T^fi, 
losing mass at a rate M, with dust expansion velocity t^oxp 
and dust-to-gas ratio the o ptical depth at wav elength A 
is approximately (see eq. 12 of lPiovan et al.ir2003l ): 



*Mka 1 



(A3) 



where nx is the extinction coefficient per unit mass and rin 
is the dusty shell internal radius. 

If the shell is optically thick to radiation of wavelength 
A, we can write 



and 



L 



1/2 



(A4) 



(A5) 



where R is the stellar radius, a is the Stefan-Boltzmann 
constant, and Td is the dust condensation temperature at 
nn. Adopting Td = 1000 K when C/O < 0.97 and Td = 
1500 K when C/O ^ 0.97, 



nn = 2.37 X 10' 



1/2 



cm (for O — rich stars) 



1/2 



Tin ~ 1.05 X 10^^ ( ) cm (for carbon stars); 
m any case (see eq. 4 of lMarigo et~alll2008l ): 



A2.1 Wind expansion velocity. 



(A6) 



(A7) 



Based on the s olutio n to the dus t- wind problem by 
lElitzur fc Ivezid (j200lh . iMarigo et aD lj2008h express the 
wind expansion velocity Hoxp as a function of M, L, ^ , and 
other dust parameters as follows: 



M 



4/3 



-1/2 



km s 



(A8) 



A = 3.08 X 10" 



B = 



rrn4 ^ 2 -1 

JdS Q*<^22 XO 

4/3 



2.28- 



Td3 



(A9) 
(AlO) 



Td3 is the dust condensation temperature in units of 10 
K; is the Planck averagj^ of the efficiency coefficient for 
radiation pressure, evaluated at Tcb, and Qv is the efficiency 
coefficient for absorption at the visual range, xo is defined 
as: 



XO 



(A12) 



with Qp{T) the Planck average of the absorption efhciency 
at temperature T, whereas 0-22 is o-gas, the dust cross-section 
per gas particle, in units of lO"^'^ cm'^, and 



' a 

= TT I - 

\2 



"■dust 2 

cm , 



(A13) 



with a the mean size of the dust grains in cm (assumed to 
be the same for all grains), and ridust and rigas the number 
densities in cm~^ of, respectively, the dust and gas particles. 

CTgas can be written as a function of <&, the dust-to-gas 
ratio, as follows: 



Pdust 
Pgas 



Pgrain^dust 



(A14) 



where pdust and pgas are the density of matter in the form of 
dust and gas, respectively, in g cm~^. j4ga3 — 4/(4XH-|-XHe) 
is the mean molecular weight of the gas in units of the H 
atom mass, mn ~ 1.674 x 10"^* g, assuming that all the 
gas is composed by H and He, with respective mass fraction 
abundances Ah and Ahc. Substituting ea. lA14l into ea. lA13l 



3 Agasmn 

2 O Pgrain 



(A15) 



The values of Qv for silicates, silicon carbide, and amor- 
phous carbon (see below. Subsection I A2 . 2 1 an d Section [A3|l 
are taken from Table 3 in IMarigo et al. I (I2OO8I ): the values of 
Q*(Teff), Qp{Tcs), and Qp{T^) (each one for all the same 
materials) are interpolated from the quantities tabulated 
there, taking temperature as the independent variable 

Clearly, eq. IA8I is valid as long as there is enough dust 
condensation to drive a n outflow. Th is condition sets the 
minimum mass- loss rate l|Elitzur fc Ivezic 2001 ): 



= 3 X lO-" 



1/2 



M, 



'"22-^'4Tfc3 



yr 



(A16) 



T he Planck mean of a function Q(A) is given by llBIanco et al.l 
Il983t) : 



< Q{T) >= 



J B(A,T)Q(A)dA 


J B(A,r)dA 



(All) 



where M_6 is the mass-loss rate in units of 10~ yr~ 
and 1/4 is the stellar luminosity in units of lO^L©. 
The parameters A and B are defined as: 



where B{\,T) is the Planck function at a temperature T. 

13 As seen in Sectional] Td = 1000 K for C/O < 0.97, and Td 

1500 K for C/O ^ 0.97. 
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Tit3 is the kinetic temperature (in units of 10'^ K) at the inner 
boundary of the dust shell, which is assumed to be equal to 
the dust condensation temperature, Tds- At smaller values 
of M, dust cannot drive a wind but may still f orm in the 
circumstellar envelope. Following iMarigo et al.l l|2008l ). we 
handle this situation by setting Uvexp = Vexp{Mmin), while 
using the actual M to calculate the opacity of the envelope. 



A2.2 Dust-to-gas ratio. 

If dust and gas share the same outflow velocity, the dust- 
to-gas ratio can also be written as the quotient between 
mass-loss in the form of dust and mass-loss in the form of 
gas: 



M - Mdus 



(A17) 



W e calculate Mdust as described by , 
(120081') ■ who follow the formalism by iFerrarottil l|2003l) 
iFerrarotti &: Gail l|2006l ). Summarizing, 



(A18) 



with the summation over several dust species, written as 



I 1 Jdust,i- 

V ^sccd / 



(A19) 



Here, Adust, i is the mean molecular weight of the dust 
species, Xgeed is the total mass fraction of the seed element 
in the circumstellar shell, and Ascad is its atomic weight. The 
seed element is the least abundant one among those needed 
to form the considered dust species, and hence limits its 
supply. Finally, /dust,i is the fraction of the seed element 
condensed into dust grains, that for each dust species is cal- 
culate d as a func t ion of M and C/O ratio using the analytic 
fits in iFerrarottil (|2003l '). as shown below. 

Stars with C/O < 0.97. In these conditions, there are 
enough oxygen atoms to form silicate- type dust, and: 



dMd 



dt 



dt 



MXsi^/sii, 



(A20) 



where Asi is the atomic weight of silicon, Asii is the effective 
molecular weight of the silicate dust, and f^n is the condensa- 
tion degree, also of the silicate dust. Given the dust mixture 
we use here for O-rich stars (see below. Section 1X3)) . the con- 
densation degree includes the contributions from amorphous 
and crystalline silicates: 



^sil — ,/warm /"cold ^~ fens /"fors 



(A21) 



where subscripts stand for warm silicate, cold silicate, en- 
statite, and forsterite dust. On the other hand, the effective 
molecular weight of the silicate dust mixture is: 



Here, Y = X/A is abundance in m olecules g ^, and lc,i = 
Yo - 2Fsi (|Ferrarotti fc Gaillbood ). Once f^u is determined 
for each star, the relative degrees of condensation of crys- 
talline and amorphous silicates are set according to the opti- 
cal depth of the envelope, as described below (Section lA3|l F°l 

Stars with C/O ^ 0.97. These are C-rich stars, since 
the dust mixture in their envelopes is dominated by car- 
bon. As stated below in Section IA3I we consider two dust 
constituents, SiC and amorphous carbon (AMC), such that: 



dMdust dA/sic , dMAMC 



dt 



dt 



dt 



(A24) 



where the first addend corresponds to the silicon carbide 
dust and the second, to the amorphous carbon dust. Fur- 
thermore, 



Marigo et all dMAMC ^ MXcfAmc 



^ii = (/. 



warm ^warm 



+ /coldj4cold + ,fonsAcns + ,/forsAfors)/,/sil.(A22) 



The degree of condensat ion of silicate ty pe dust is 
fo und from the analyt ic fit in iFerrarott J (12003), used also 
bylMarigoelaD (|2008l ): 



/sii = 0.8 



Af_6 + 5 



Yc,i - Yc 



(A23) 



dt 



(A25) 



Xc is the carbon abundances by mass; we take the 
current atmospheric value, solar scaled, or Xc = 
{C/0)Xo,q{Z/Zq), where Xq.q is th e solar mass frac- 
tion a bundance of oxygen taken from iGrevesse fc Sauval 

Finally, the deg ree of condensation of carbon is found 
from IFerrarottil ([2003)13 



/amc ~ 0.5 



M- 



Af_6 + 5 



f- 

Vo 



0.97 



(A26) 



Once /amc is found, the relative degree of condensation of 
SiC dust is set according to the opacity of the envelope, as 
described below (Section |A3|I . 



A2.3 Extinction coefficient. 

In order to determ ine ka, we follow the procedure by 
IPiovan et ahi l|2003l ). The mass extinction coefficient can be 
written as: 



^ ^dust,z'^dust,^ 
pdust 



(A27) 



with the summation over all types of grains in the dust mix- 
ture, and where the ith type of grain has number density 
Wdust.i cm~'^, and cross-section for radiation-dust interac- 
tions (Tdust,i cm^; Pdust (g cm^) is, again, the density of mat- 
ter in the form of dust, and (Tdust,i(a) ~ i^a^Qa^t{i), where 
Qc^tii) are the extinction coefflcients. If mdust, i is the mass 
of the ith type of dust grains, and we introduce the mass 
abundance x^ = ndust,imdust,i/Pdust: 



E 



?^dust,i 



(A28) 



' For the cold and warm amor phous silicates we us e 
the following formulae, respectively llDorschner et al.l Il995h : 

Mgo.8Fei.2Si04; Mgn.4F eo.6Si03. 

For consistency with iMarigo fc Girardil l|20o3), and with the 
different works by Ferrarotti and collaborators. 
22 Ferrarotti writes (C/O - 1), but we change it to (C/O - 0.97), 
to avoid the possibility of negative degrees of condensation. 
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Going back to eq. IA3I one can see that dust optical 
depth is a function of ka, but then is a function of the 
optical depth, because the dust composition is itself a func- 
tion of T\. Consequently, an iterative procedure is needed to 
determine the dust optical depth. 



A3 Implementation of the envelopes with 
DUSTY. 

The DUSTY code includes the dependence of the wind speed 
on stellar luminosity and metallicity (through the gas-to- 
dust ratio of the envelope), and the drift speed between the 
dust and the gas. The code outputs spectral shapes, that can 
then be scaled to energy flux by multiplying by the stellar 
luminosity, and dividing by ATvd^ . DUSTY offers a broad 
choice of input parameters, specifically in regard to the dust 

chemical composition and grain size distribution. 

We adopt he re the dust mixtures introduced by ISuhl 
1 19991 . 12OOOI . 120021 ') ■ with the aim of fitting o bservations of 
individual stellar spectra, and later used by IPiovan et al.l 
(|2003l l in SSP models. Suh has calculated optical constants 
for siliceous and carbonaceous compounds that are consis- 
tent with the Kramers-Kronig dispersion relations and, at 
the same time, yield models of dusty envelopes that fit ob - 
served properties of AGB stars. For O-rich stars. [Suhl (|2002l ') 
proposes combinations of both amorphous and crystalline 
silicate grains (enstatite, MgSiOs, and forsterite, Mg2Si04), 
in proportions that depend on the optical depth of the dusty 
shell: (a) for stars with the lowest mass-loss rates and thin 
shells, where the 10 ^m silicate feature is observed in emis- 
sion (no 3), a blend of 90% warm amorphous silicate, 
5% enstatite, and 5% forsterite grains; (6) for stars with 
low mass-loss rates and moderately optically thick shells (3 

< Tio ^ 15), 90% cold amorphous silicate, 5% enstatite, 
and 5% forsterite grains; (c) for stars with high mass-loss 
rates and optically thick shells (no > 15), 80% cold amor- 
phous silicate, 10% enstatite, and 10% forsterite. Dust opac- 
ity functions for these 3 cases a re sh o wn in Fig. lAll 

In the case of C-rich stars, 'Suh' (2000) prescribes mix- 
tures of amorphous carbon (AMC) and silicon carbide (SiC), 
again in ratios that depend on the optical depth of the dusty 
shell, (a) For C-rich stars with optically thin shells (no ^ 
0.15), 20% of SiC grains are needed to fit the strong 11 
^m feature in observed spectra; (b) 10% of SiC grains are 
required for shells with intermediate optical depths (0.15 

< no ^ 0.8), and (c) no SiC grains are necessitated at larger 
optical depths, for which the 11 /im feature is either weak or 
absent. Dust opacity functions for the 3 cases encountered 
for C-rich stars are displayed in Fig. IA2I 

Optica l constants for the differ ent g r ains a re taken from 

ISuhl (|l999l, amorpho u s silic ates), ISuhl l|2000l , amorphous 
carbo n) , Jaeger et al.l l| 19981 . enstatite) , and iFabian et al.l 
l|200ll . forsterite -actua lly, olivin e : MR i.9Feo.iSi04). The 
data for a-SiC come from lPegourTi l|l98d ). and are included 
in DUSTY's built-in library of optical constants. 

We take a uniform dust grain size of 0.1 /im, follow- 
ing [Piwa^^er^l] (|2003l ). In general, for stars with M 7^ 0, 
we use DUSTY's numerical solution for radiatively driven 
winds, extending to a distance 10* times the inner radius. 
For T ^ 80, however, and due to numerical difficulties of 
the program, we assume a shell density distribution that 
falls off as r~^. As for the incident radiation (that DUSTY 



surmises comes from a point source at the centre of the den- 
sity distribution), we use the same stellar SEDs as the BC03 
and CB09 standard models (adapted according to stellar pa- 
rameters for modified mass-loss-rates) , and for the dust tem- 
perature on the envelope inner boundary we input 1000 K 
for O -rich stars, and 1500 K fo r carbon stars (iMarigo et al.l 
2008hFI For consis tency with [Su3 (|l999l . |2000| . l2002h "and 
Piovan et ahl (120031 ) . we evaluate k,\ and ta at lO^m, and 
naturally all our calculations with DUSTY take the same A 
as their reference wavelength. 



REFERENCES 

Ajhar, E. A. & Tonry, J. L. 1994, ApJ, 429, 557 

Ajhar, E. A., Lauer, T. R., Tonry, J. L., Blakeslee, J. P., 

Dressier, A., Holtzman, J. A., & Postman, M. 1997, AJ, 

114, 626 

Alongi, M., Bertelli, G., Bressan, A., Chiosi, C, Fagotto, 
F., Greggio, L., & Nasi, E. 1993, A& AS., 97, 851 

Asplund, M., Grevesse, N., & Sauval, A. J. 2005, Cosmic 
Abundances as Records of Stellar Evolution and Nucle- 
osynthesis, 336, 25 

Baud, B., & Habing, H. J. 1983, A& A, 127, 73 

Bertelh, G., Girardi, L., Marigo, P., & Nasi, E. 2008, A& 
A, 484, 815 

Blakeslee, J. P., et al. 2009, ApJ, 694, 556 

Blanco, A., Falcicchia, G., & Merico, F. 1983, Ap&SS, 89, 
163 

Bloecker, T. 1995, A& A, 297, 727 

Bowen, G. H., & Willson, L. A. 1991, ApJL, 375, L53 

Bressan, A., Fagotto, F., BerteUi, G., & Chiosi, C. 1993, 

A& AS., 100, 647 
Bressan, A., Granato, G. L., & Silva, L. 1998, A& A, 332, 

135 

Bruzual A., G. 2002, in Extragalactic Star Clusters, lAU 
Symposium Ser., Vol. 207, eds. D. Geisler, E. K. Grebel, 
& D. Minniti (Provo:ASP), 616 

Bruzual, G., 2007, in Proceedings of the lAU Symposium 
No. 241 "Stellar populations as building blocks of galax- 
ies", eds. A. Vazdekis and R. Peletier, Cambridge: Cam- 
bridge University Press, 125 (arXiv:astro-ph 0703052) 

Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 

Chariot, S., & Bruzual, G. 2009, in preparation 

Cantiello, M., Blakeslee, J. P., Raimondo, G., Mei, S., Bro- 
cato, E., & Capaccioh, M. 2005, ApJ, 634, 239 

Cantiello, M., Raimondo, G., Brocato, E., & Capaccioli, 
M. 2003, AJ, 125, 2783 

Cantiello, M., Blakeslee, J., Raimondo, G., Brocato, E., & 
Capaccioh, M. 2007, ApJ, 668, 130 



For the most massive stars (M > 4Mq), both C-rich and 
O-rich, occasionally there are extremely short (10 — 100 yr long) 
phases with no mass-loss that follow mass-losing stages with a du- 
ration of ~ 10® yr. We worried that standard (i.e., Schultheis et 
al. or Westera et al., as appropriate) spectra would not provide a 
realistic modelling of such phases, and performed the experiment 
of assigning to them, instead, the spectra of the immediately pre- 
ceding stage. The results of both procedures (standard spectra 
and preceding stage spectra) were indistinguishable for all prac- 
tical purposes. 



Mass-loss tracers 29 




Figure Al. Absorption efficiency functions of different mixtures of siliceous dust grains, for envelopes of O-rich stars with Vcirious optical 
depths. Dotted-line: no ^ 3; solid-line: 3 < no ^ 15; dashed-line: no > 15 (see text). 



Cantiello, M., Raimondo, G., Blakeslee, J. P., Brocato, E., 

& Capaccioli, M. 2007, ApJ, 662, 940 
Corvino, Valls-Gabaud, Luridiana, & Mas-Hesse, J. M. 

2002, A&A, 381, 51 
Cervifio, M., & Luridiana, V. 2006, A& A, 451, 475 
Chabrier, G. 2003, PASP, 115, 763 

Cioni, M.-R., et al. 2000, VizieR Online Data Catalog, 

2228, 

Cohen, J. G. 1982, ApJ, 258, 143 

de Jager, C, Nieuwenhuijzen, H., & van der Hucht, K. A. 

1988, A& AS., 72, 259 
Demleitner, M., Accomazzi, A., Eichhorn, G., Grant, C. S., 

Kurtz, M. J., & Murray, S. S. 2001, ASP Conf. Ser. 238: 

Astronomical Data Analysis Software and Systems X, 238, 

32 

Dorschner, J., Begemann, B., Henning, T., Jaeger, C, & 

Mutschke, H. 1995, A& A, 300, 503 
Elitzur, M., & Ivezic, Z. 2001, MNRAS, 327, 403 
Elson, R. A. W. & Fall, S. M. 1985, ApJ, 299, 211 

1988, AJ, 96, 1383 

Epchtein, N., et al. 1997, The Messenger, 87, 27 

Fabian, D., Henning, T., Jagor, C, Mutschke, H., 

Dorschner, J., & Wehrhan, O. 2001, A& A, 378, 228 
Fagotto, F., Bressan, A., Bertelli, G., & Chiosi, C. 1994, 

A& AS., 104, 365 



1994, A& AS., 105, 29 

Ferrarese, L. et al. 2000, ApJS, 128, 431 
Ferraro, F. R., Fusi Pccci, F., Testa, V., Groggio, L., Corsi, 
C. E., Buonanno, R., Torndrup, D. M., & Zinncckcr, H. 

1995, MNRAS, 272, 391 

Ferrarotti, A. S., & Gail, H.-P. 2001, A& A, 371, 133 
Ferrarotti, A. S. 2003, Ph.D. Thesis 
Ferrarotti, A. S., & Gail, H.-P. 2006, A& A, 447, 553 
Frogel, J.A., Mould, J., & Blanco, V.M., 1990, ApJ, 352, 
96 

Gail, H.-R, & Sodlmayr, E. 1986, A& A, 161, 201 
Gallart, C, et al. 2003, AJ, 125, 742 

Girardi, L., Chiosi, C, Bertelli, G., & Bressan, A. 1995, 

A&: A, 298, 87 

Girardi, L., Bressan, A., Chiosi, C, Bertelli, G., & Nasi, E. 

1996, A&: AS., 117, 113 

Girardi, L., & Marigo, P. 2007, A& A, 462, 237 
Gonzalez, R. A., Liu, M. C, & Bruzual A., G. 2004, ApJ, 
611, 270 

2005, ApJ, 621, 557 

Gonzalez- Lopezlira, R. A., Albarran, M. Y., Mouhcine, M., 

Liu, M. C, Bruzual-A., G., & de Batz, B. 2005, MNRAS, 

363, 1279 

Gonzalez-Lopezlira, R. A., & Buzzoni, A. 2007, in Ul- 
traviolet properties of evolved stellar populations, eds. 



30 Gonzalez- Lopezlira et al. 




Figure A2. Absorption efficiency functions of different mixtures of carbonaceous dust grains, for envelopes of C-ricli stars with various 
optical depths. Dotted-line: rio ^ 0.15; solid-line: 0.15 < rio ^ 0.8; dashed-line: tiq > 0.8 (see text). 



M. Chavez, E. Bertone, D. Rosa-Gonzalez, & L. H. 
Rodriguez-Merino, Astrophysics and Space Science Pro- 
ceedings (Dordrecht:Springer), in press (arXiv:astro-ph 
0708.0193) 

Grevesse, N., & Sauval, A. J. 1998, Space Science Reviews, 
85, 161 

Groenewegen, M. A. T., & de Jong, T. 1994, A& A, 288, 
782 

Groenewegen, M. A. T., Smith, C. H., Wood, P. R., Omont, 

A., & Fujiyoshi, T. 1995, ApJL, 449, L119 
Iben, I., Jr. 1984, ApJ, 277, 333 
Iben, I., Jr., & Renzini, A. 1983, ARA&A, 21, 271 
Ivezic, Z., & Elitzur, M. 1997, MNRAS, 287, 799 
Ivezic, Z., Nenkova, M., & Elitzur, M. 1999, User Man- 
ual for DUSTY, University of Kentucky Internal Report 
( |http: //www.pa.uky ■ edvi/~moshe/dusty[ ) 
Jaeger, C., Molster, F. J., Dorschner, J., Henning, T., 
Mutschke, H., & Waters, L. B. F. M. 1998, A& A, 339, 
904 

Jensen, J. B., Tonry, J. L., Barris, B. J., Thompson, R. I., 
Liu, M. C., Rieke, M. J., Ajhar, E. A., & Blakeslee, J. P. 
2003, ApJ, 583, 712 

Kalirai, J. S., Hansen, B. M. S., Kelson, D. D., Reitzel, 
D. B., Rich, R. M., & Richer, H. B. 2008, ApJ, 676, 594 

Kalirai, J. S., Richer, H. B., Reitzel, D., Hansen, B. M. S., 



Rich, R. M., Fahlman, G. G., Gibson, B. K., & von Hippel, 
T. 2005, ApJL, 618, L123 
Knapp, G. R. 1985, ApJ, 293, 273 

Kuntschner, H. 1998, Ph.D. thesis. University of Durham, 
UK. 

LauQon, A., & Mouhcine, M. 2002, A& A, 393, 167 
Lebzeher, T., Posch, T., Hinkle, K., Wood, P. R., & Bouw- 

man, J. 2006, ApJL, 653, L145 
Lepine, J. R. D., Ortiz, R., & Epchtein, N. 1995, A& A, 

299, 453 

Liu, M. C., Chariot, S., & Graham, J. R. 2000, ApJ, 543, 
644 

Liu, M. C, Graham, J. R., & Chariot, S. 2002, ApJ, 564, 
216 

Lyubimkov, L. S., Rostopchin, S. I., Rachkovskaya, T. M., 

Poklad, D. B., & Lambert, D. L. 2005, MNRAS, 358, 193 
Maraston, C. 1998, MNRAS, 300, 872 

2005, MNRAS, 362, 799 

Maraston C, Daddi E., Renzini A., Cimatti A., Dickinson 

M., Papovich C., Pasquah A., & Pirzkal, N. 2006, ApJ, 

652, 85; (Erratum: 2007, ApJ, 656, 1241) 
Marigo, P., Bressan, A., & Chiosi, C. 1996, A& A, 313, 545 
Marigo, P., & Girardi, L. 2007, A& A, 469, 239 
Marigo, P., Girardi, L., Bressan, A., Groenewegen, 

M. A. T., Silva, L., & Granato, G. L. 2008, A& A, 482, 



Mass-loss tracers 31 

883 

Marin- Franch, A., & Aparicio, A. 2006, A& A, 450, 979 
Mei, S., Quinn, P. J., & Silva, D. R. 2001, Ak A, 371, 779 
Mei, S., et al. 2005, ApJS, 156, 113 

2005, ApJ, 625, 121 

2007, ApJ, 655, 144 

Mighell, K. J., Rich, R. M., Shara, M., & Fall, S. M. 1996, 

AJ, 111, 2314 
Mokiem, M. R., et al. 2007, A& A, 473, 603 
Mouhcinc, M., Gonzalez, R. A., & Liu, M. C. 2005, MN- 

RAS, 362, 1208 
Mouhcine, M., & Langon, A. 2003, A& A, 402, 425 
Pahrc, M. A., & Mould, J. R. 1994, ApJ, 433, 567 
Pegouric, B. 1988, A& A, 194, 335 

Piovan, L., Tantalo, R., & Chiosi, C. 2003, A& A, 408, 559 
Raimondo, G. 2009, ApJ, 700, 1247 

Raimondo, G., Brocato, E., Cantiello, M., & Capaccioli, 

M. 2005, AJ, 130, 2625 
Reimers, D. 1975, Mcmoires of the Societe Royale des Sci- 
ences de Liege, 8, 369 

1977, A& A, 61, 217 

Renzini, A., & Buzzoni, A. 1986, Spectral Evolution of 

Galaxies, 122, 195 
Renzini, A., & Voli, M. 1981, A& A, 94, 175 
Russell, S. C., & BossoU, M. S. 1989, ApJS, 70, 865 
Santos, J. F. C., & Frogcl, J. A. 1997, ApJ, 479, 764 
Schultheis, M., Aringcr, B., Jorgonscn, U. G. Lebzelter, 

T., & Ploz, B. 1999, Abstract of the 2nd Austrian ISO 

workshop: Atmospheres of M, S, and C Giants, eds. J. 

Hron & S. Hofner, 93 
Searle, L., Wilkinson, A., & Bagnuolo, W. G. 1980, ApJ, 

239, 803 

Skrutskie, M. F. et al. 1997, ASSL Vol. 210: The Impact of 

Large Scale Near-IR Sky Surveys, 25 
Steflen, M., Szczerba, R., Men'shchikov, A., & Schoen- 

berner, D. 1997, A& AS., 126, 39 
Suh, K.-W. 2002, MNRAS, 332, 513 

2000, MNRAS, 315, 740 

1999, MNRAS, 304, 389 

Srinivasan, S., et al. 2009, AJ, 137, 4810 

Tonry, J. L., Ajhar, E. A., & Luppino, G. A. 1990, AJ, 100, 

1416 

Tonry, J. L., Blakeslee, J. P., Ajhar, E. A., & Dressier, A. 

1997, ApJ, 475, 399 
Tonry, J. L., Dressier, A., Blakeslee, J. P., Ajhar, E. A., 

Fletcher, A. B., Luppino, G. A., Metzger, M. R., & Moore, 

C. B. 2001, ApJ, 546, 681 
Tonry, J. & Schneider, D. P. 1988, AJ, 96, 807 
van den Bergh, S. 1981, A&AS, 46, 79 
van Loon, J. T. 2000, A& A, 354, 125 

2006, Stellar Evolution at Low Metalhc- 

ity: Mass Loss, Explosions, Cosmology, 353, 211 
van Loon, J. T., Cioni, M.-R. L., Zijlstra, A. A., & Loup, 

C. 2005, A& A, 438, 273 
Vassiliadis, E. & Wood, P. R. 1993, ApJ, 413, 641 
Wachter, A., Schroder, K.-P., Winters, J. M., Arndt, T. U., 

& Sedlmayr, E. 2002, A& A, 384, 452 
Westera, P., Lejeune, T., Buser, R., Cuisinier, F., & 

Bruzual A., G. 2002, A& A, 381, 524 
Willson, L. A. 2000, ARA&A, 38, 573 
Worthey, G. 1993a, ApJ, 415, 91 
1993b, ApJ, 418, 947 



